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PREFACE 

This report presents the results of studies conducted during the 
period July 19, 1969 — July 19, 1970, under NASA Research Contract 
NAS 8-21432, "Lunar Surface Engineering Properties Experiment Definition." 
This study was sponsored by the Lunar Exploration Office, NASA Head- 
quarters, and was under the technical cognizance of Dr. N. C. Costes, 

Space Science Laboratory, George C. Marshall Space Flight Center. 

The report reflects the combined effort of five faculty investiga- 
tors, a research engineer, a project manager, and eight graduate research 
assistants, representing several engineering and scientific disciplines 
P ert i nen t to the study of lunar surface material properties. James K. 
Mitchell, Professor of Civil Engineering, served as Principal Investigator 
and was responsible for those phases of the work concerned with problems 
relating to the engineering properties of lunar soils and lunar soil 
mechanics. Co-investigators were William N. Houston, Assistant Professor 
of Civil Engineering, who was concerned with problems relating to the 
engineering properties of lunar soils; Richard E. Goodman, Associate 
Professor of Geological Engineering, who was concerned with the engineer- 
ing geology and rock mechanics aspects of the lunar surface; Paul A. 
Witherspoon, Professor of Geological Engineering, who was concerned with 
fluid conductivity of lunar surface materials in general; Franklin C. 
Hurlbut, Professor of Aeronautical Science, who was concerned with 
experimental studies on fluid conductivity of lunar surface materials; 
and D. Roger Willis, Associate Professor of Aeronautical Science, who 
conducted theoretical studies on fluid conductivity of lunar surface 
materials. Dr. Karel Drozd, Assistant Research Engineer, performed 
laboratory tests and analyses pertinent to the development of a borehole 
jack for determination of the in situ characteristics of lunar soils 
and rocks; he also helped in the design of the borehole jack. H. Turan 
Durgunoglu, H. John Hovland, Laith I. Namiq, Parabaronen Raghuraman, 

James B. Thompson, Donald D. Treadwell, C. Robert Jih, Suphon Chirapuntu, 
and Tran K. Van served as Graduate Research Assistants and carried 
out many of the studies leading to the results presented in this 
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report. Ted S. Vinson, Research Engineer, served as project manager 
until May 1970, and contributed to studies concerned with lunar soil 
stabilization. H. John Hovland served as project manager after May 
1970, and contributed to studies concerned with soil property evaluation 
from lunar boulder tracks. 


Ultimate objectives of this project were: 


1) Assessment of lunar soil and rock property data using information 
obtained from Lunar Orbiter, Surveyor, and Apollo missions. 

2) Recommendation of both simple and sophisticated in situ testing 
techniques that would allow determination of engineering 
properties of lunar surface materials. 

3) Determination of the influence of variations in lunar surface 
conditions on the performance parameters of a lunar roving 
vehicle. 


4) Development of simple means for determining the fluid 
conductivity properties of lunar surface materials. 


5) 


Development of stabilization techniques for use in loose, 
unconsolidated lunar surface materials to improve the 


j-'cj- -L wi.nianvc o f ouK^ii luat^rials in luMuiT engineering application. 


The scope of specific studies conducted in satisfaction of these objectives 
is indicated by the following list of contents from the Detailed Final 
Report which is presented in four volumes. The names of the investigators 
associated with each phase of the work are indicated. 


VOLUME I 


MECHANICS, PROPERTIES, AND STABILIZATION OF LUNAR SOILS 

1. Lunar Soil Simulant studies 

W. N . Houston, L. I. Namiq, J. K. Mitchell, and D. D. Treadwell 

2. Determination of In Situ Soil Properties Utilizing an Impact 
Penetrometer 

J. B. Thompson and J. K. Mitchell 

3. Lunar Soil Stabilization Using Urethane Foamed Plastics 
T. S. Vinson, T. Durgunoglu, and J. K . Mitchell 

4. Feasibility Study of Admixture Soil Stabilization with Phenolic 
Resins 

T . Durgunoglu and J. K. Mitchell 


VOLUME II 


MECHANICS OF ROLLING SPHERE-SOIL SLOPE INTERACTION 
H. J. Hovland and J. K. Mitchell 


1 . Introduction 

2. Analysis of Lunar Boulder Tracks 

3. Model Studies of the Failure Mechanism Associated with a Sphere 
Rolling Down a Soil Slope 

4. Pressure Distribution and Soil Failure Beneath a Spherical Wheel 
in Air-Dry Sand 

5. Theoretical Studies 

6. Rolling Sphere Experiments and Comparison with Theoretical 
Predictions 

7 . Utilization of Developed Theory 

8. Conclusions and Recommendations 
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1. Summary of Previous Work 

R. E. Goodman, T. K. Van, and K. Drozd 

2. An Experimental Study of the Mechanism of Failure of Rocks 
Under Borehole Jack Loading 

T. K. Van and R. E. Goodman 

3. A Borehole Jack for Deformability , Strength, and Stress 
Measurements in a 2-inch Borehole 

R. E. Goodman, H. J . Hovland, and S. Chirapuntu 
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FLUID CONDUCTIVITY OF LUNAR SURFACE MATERIALS 

1. Studies on Fluid Conductivity of Lunar Surface Materials — 
Theoretical Studies 

P. Raghuraman and D. R. Willis 

2. Studies on Fluid Conductivity of Lunar Surface Materials — 
Experimental Studies 

F. C. Hurlbut, C. R. Jih, and P. A. Witherspoon 
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VOLUME II 


Mechanics of Rolling Sphere — Soil Slope Interaction 
Chapter 1. INTRODUCTION 

This report presents the results of studies related to rolling- 
sphere (rolling boulder) soil-slope interaction. It describes all work 
done on this phenomenon during the last 2 years at Space Sciences 
Laboratory, University of California. This work was also presented as 
a Ph.D. dissertation by H. John Hovland to the University of California. 

The motivation for investigating the relationship between soil 
properties and the track left by a rolling sphere was provided by the 
observation of boulder tracks on the surface of the moon in Lunar Orbiter 
photographs. A typical lunar boulder track is shown in Figure 1-1. 

It was soon recognized that a relationship must exist between the 
size and the track of the boulder, and that this relationship would 
include both soil and boulder properties. Hence, some lunar boulders 
were studied (Filice, 1967; Nordmeyer and Mason, 1967; Moore and Martin, 
1967; Eggleston et al., 1968) in an effort to determine the static 
bearing capacity of lunar surface soil. We have been investigating in 
more detail the possibilities for deducing soil strength parameters from 
the boulder tracks. 

A summary of the work done by our group on the study of lunar boulder 
tracks during the year 1967-68 was presented in the final report for 
Contract NSR 05-003-189 (Mitchell et al., 1969). In that report, several 
methods for analyzing boulder-track data were considered, each subject 
to limitations and giving somewhat different results. It was recommended 
that boulder-track phenomena be further studied, and it was noted also 
that if variability of the lunar surface is to be determined, it is 
important to use the same method of analysis throughout. 

Studies during 1968-69 consisted of reviewing Orbiter photography 
for locating suitable boulder tracks for analysis, study of site geology 
for selected boulder tracks, and static analysis of the boulder tracks 







using bearing capacity theory* These investigations were reported in 
detail by Hovland and Mitchell (1970) . The results suggested a vari- 
ability of lunar soils. Possible causes of such variability were 
considered. It was also recognized that the theory used neglected the 
dynamic aspect of the problem. It was recommended that an improved 
theory or method be developed. The pertinent results and conclusions 
from this work are presented in Chapter 2. 

The research for this report has consisted of an in-depth investi- 
gation of the fundamental nature of the rolling sphere-soil slope 
interaction phenomenon. Specific objectives of this research were: 

1. To develop an understanding of the soil deformation mode under 
the action of a rolling sphere. 

2. To develop an improved theory or method, based on the actual 

soil failure mechanism, which would provide a remote recon- 

# t > 
naissance technique for study of soil conditions using 

boulder track observations. 

The rolling sphere-soil slope interaction phenomenon was studied 
systematically to meet the stated objectives. The failure mechanism 
investigated using models and by testing an instrumented spherical wheel. 
The model studies are reported in Chapter 3. The instrumented spherical 
wheel was specifically designed to measure contact pressure, but it also 
provided much information on the failure mechanism, as reported in Chapter 
4. Based on an understanding of the failure mechanism, theory was modi- 
fied and developed as reported in Chapter 5. Further tests were conducted 
to study the dynamic aspects of the problem and the combined effect of 
several parameters. These tests included rolling some two hundred spheres 
down sand slopes. Films were taken of the rolling spheres, and the tracks 
were measured. This experimental work as reported in Chapter 6 provides 
the basis for comparison with theory. 

Implications of the results and reevaluation of the lunar boulder 
tracks analyzed using the proposed method are discussed in Chapter 7. 
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Chapter 2. ANALYSIS OF LUNAR BOULDER TRACKS 

Among the conspicuous and interesting features on the lunar surface 
are large boulders and the tracks they left as they rolled down slopes. 

These features were observed early on photographs provided by lunar 
orbiters. The tracks appear to be of three different types; i.e., 

(1) continuous tracks suggestive of a spherical boulder rolling uniformly 
down the slope, (2) segmented tracks suggestive of bouncing, and (3) 
relatively short tracks suggestive of plowing. 

At the outset of this investigation sixty-nine boulder tracks from 
nineteen different locations on the moon, as shown in Figure 2-1, were 
analyzed. These locations include upland, maria, and perhaps intermediate 
terrain. Only tracks that appeared continuous (implying a relatively 
spherical boulder and uniform rolling) were selected. 

A complete report including consideration of the geology of boulder 
track locations and detailed discussion of the implications of the results 
was presented previously (Hovland and Mitchell, 1970). A summary of this 
study is given here. 

METHOD OF ANALYSIS 

Geometrical Relations of Sphere and Track 

A boulder rolling down a slope where the soil fails in general shear 
would leave a track with a raised rim, (Figure 2-2). Such rims have 
been observed on many lunar boulder tracks. For the purpose of the 
present analysis , preliminary theory can be developed for a somewhat more 
idealized situation, assuming a sphere -track geometry as shown on Figure 2-3. 

From Figure 2-3 it may be seen that the track depth will be given by 
z = r(l - cos0) * r (l - cos £sin * jj' J ) , (2-1) 


where 


D = 2r, or the sphere diameter; and w is the track width 




FRAME LOCATIONS 
A ORBITER II 
□ ORBITER III 
O ORBITER V 


Fig. 2-1. Locations of boulder tracks analyzed. 




FAILURE SURFACE BREAKOUT 


TOP VIEW 



FRONT VIEW 



SIDE VIEW 


Fig. 2-2. Boulder rolling down slope. 



SIDE VIEW 



TOP VIEW (NORMAL TO SLOPE) 



Fig. 2-3. Geometrical relations of sphere and track. 
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The semicircular soil-sphere contact area may be represented by an 
equivalent rectangular area defined by 


giving 


2b 1 = S- I 

4 2 


f 



0.44 w . 


( 2 - 2 ) 

(2-3) 


Modified Bearing Capacity Theory 

A general bearing capacity equation for a strip footing (Leonards, 1962) 
is 

V 

q = — N y + cN c + q'N q • (2-4) 

This equation assumes a homogeneous-isotropic soil; other assumptions are 
discussed in Chapter 5, pages 95 and 101. For a rectangular footing this 
equation may be modified to 

V 

q * — Vy * c Vc + q ' Vq • <2 " 5) 

In these equations 

q = unit bearing capacity 

Y = soil density in earth gravity 
s 

b = width of footing 
c = soil cohesion 
q' = surcharge 

s , s , s = shape factors, and 
y c q 

N , N , N = bearing capacity factors which depend on the soil 
Y c q 

friction angle, <j>. 

Skempton (1951) indicated that, for <|> = 0, the value of s can be 

c 

taken as (1 + 0.2 b/L) , where L is the length of a rectangular footing. 

For 4> > 0, the vaue of s c would probably not be significantly different. 
Meyerhof (1951) proposed that, for 4> = 30°, s^ equals approximately 



(1 + 0.2 b/L) . The friction angle for lunar soil is probably close enough 
to 30° to justify the use of this value for s . The shape factor, s . 

q ^ y 

is given by (1 - 0.3 b/L) according to Lundgren and Hansen (1955) and 
Hansen (1957) . 


Substituting these shape factors into the bearing capacity equation 
we get 



0.3 |)b y + c (l + 0.2 |)b c + q' 



(2-6) 


for the sphere b/L = 1/2 and b = 0.444 w. If, in addition, an average 
surcharge depth is taken to be z/2. 


q = °*^ 5 (0.444 w)y N + 1.1 cN + y zN . 

2 s y c 2 's q (2-7) 

A convenient equation in dimensionless form results by dividing both 

sides by (wy ) to give 
s 


-J- = 0.188 N + 1.1 N + 0.55 ( N . 

wy s y ywYg/ c Vw/ q 


(2-8) 


Defining q as the unit bearing capacity in earth gravity and Y as 
e s 

the unit weight of soil in earth gravity, the equation for the unit bearing 
capacity of a sphere in the earth gravity field becomes 


~r = 0.188 N + 1.1 + 0.55 (-) N . (2-9) 

WY S y v^s/ c 'w/ 3 

This equation is readily adapted for estimating the bearing capacity 
on the moon by noting that, for a given soil mass density, the unit weight 
on the moon will be reduced by a factor of 6. Thus, if q^ designates the 
bearing capacity on the moon. Equation (2-9) becomes 


— = 0.0314 N y + 1.1 
s 



+ 0.0916 N . 

Vw / q 


( 2 - 10 ) 
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The influence of the slope angle, a, can be incorporated by using 
Meyerhof's bearing capacity factors (Meyerhof, 1951) . For the purpose 
of this study, the charts for Meyerhof's bearing capacity factors were 
enlarged and are presented on Figures 2-4, 2-5, and 2-6. 

An expression for the average bearing pressure under a sphere can be 
developed from the ratio of sphere weight to bearing area. This gives 


4/3TTr 3 Y r 

IL iii 

2 4 


32 

3 



4 

3 


Dy 

r 

(w/D) 2 


( 2 - 11 ) 


where 

r = sphere radius 
D = sphere diameter, and 
Y r = rock or sphere density in earth gravity . 

Dividing both sides of this equation by wy we have 

S 9 



4 W^s > 
3 (w/D) 3 


( 2 - 12 ) 


Again, for the lunar gravity field, lunar soil unit weight equals v /6 

's' 

Then the average bearing pressure for a sphere in lunar gravity is given by 


2 Y r/ Y s 

^s 9 (w/D) 3 ' (2-13) 


Table 2-1 lists some possible rock and soil combinations that might 
be found on the lunar surface. 




Fig. 2-4. General bearing capacity factor N £ for strip 
foundation (after Meyerhof, 1951). 
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Table 2-1. Possible Lunar Rocks and Soils. 


Description 

Likelihood 

Rock type 

Y r 

1 gm/cc 

Soil type 

*s 

gm/cc 

Y r 

Y s 

of ratio 
occurrence 

Basaltic, vesicular 

2.7 

Mixed and maria soil 

1.6 

1.7 

Probably 






common 

Basaltic, average 

porosity (Surveyor 
results) 

3.1 

Mixed and maria soil 
(Surveyor results) 

1.5 

2.0 

Probably 

common 

Porous volcanic rock 

1.8 

Relatively dense 
basaltic soil 

1.8 

1 

Possible 

Basaltic, average 
porosity 

3.1 

Volcanic ash or 
porous material 

1.1 

3 

Possible 

Extra dense volcanic 
or meteoritic rock 

4.2 

Light volcanic ash or 
porous material* 

0.7 

6 

Probably 

rare 

Pumice-like rock 

0.7* 

Dense basaltic soil 

2.3 

0.3 

Probably 
very rare 


*Bendix Corporation 1968 


Equation (2-13) was solved using the density ratios listed in 
Table 2-1 and gave the results plotted in Figure 2-7. 

The angle of internal friction, <f>, can now be determined by using 
Equation (2-10) and either Equation (2-13) or Figure 2-7, provided c 
is known or can be estimated. It is to be noted that Equation (2-10) is 
an upper bound to bearing capacity because it gives a solution based on 
the maximum resistance available for given values of c and <{>. In prin- 
ciple, Equation (2-13) is an exact value of bearing capacity, since for 
any lower value of q m the track width would be greater, and for any 
greater value of q m the track width would be smaller. 

SELECTION OF BOULDERS AND MATERIAL PROPERTIES 

The procedure followed for study of different boulder tracks con- 
sisted of (1) locating suitable boulder tracks on high resolution Lunar 
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Fig. 2-7. 
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Graph of Equation (2-13). 





Orbiter photographs, (2) determining the scale of the photograph or 
frame, (3) measuring the boulder and track dimensions, and (4) estima- 
ting the slope angle for selected portions of the track. Some of the^ 

v 

above phases of the analysis are described in greater detail below. 

Scale of the Photographs 

The Orbiter Supporting Data* for the Lunar Orbiter missions give 
the latitude and longitude of 44 equally spaced angles along the photo 
frame periphery. Points 1, 12, 23, and 34 correspond to the corners of 
the frame, and can be easily identified. Hence, the distance from one 
corner to another can be measured on the frame, and the corresponding 
ground distance can be calculated from the latitude and longitude of 
the corners given in the supporting data. Thus, the scale to be used is 

, Distance on Frame 

Scale = 

Distance on Ground 


The distance on the ground is either given in the supporting data 
or can be determined from geometrical relationships of a sphere. The 
length of any circular arc on a sphere is given by 


AB 


TTR 

m 

180 


AB 


(2-14) 


where 


L = arc distance between points A and B 
AB 

R = radius of the moon = 1740 km 
m 

0 = angle between A and B. 

AB 


The geometrical relationships are illustrated on Figure 2-8. By the 
Pythagorean theorem, neglecting surface curvature, we have for the 
central angle 



* Revised data dated 2-5-69 give the latitude and longitude of addi- 
tional points including the corners of the frame (Boeing, 1969) . 




Fig. 2-8. Geometry for determining arc length. 
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In these equations 







= longitude of point A, 

= latitude of point A, 

= longitude of point B, 

= latitude of point B, and 


0 LO is a< ^ jus ted by cos — since a degree longitude is smaller 

near the poles than near the equator. A more correct expression for 
the distance which accounts for surface curvature (Roggeveen and 

Goodman 1968), is 


where 



9 ' 



For central angles less than 3 degrees, cos 6' & 0.999. Therefore, 
for most of the Orbiter high-resolution photographs where boulder 
tracks are found. Equation (2-15) is sufficiently accurate. 

The Orbiter Supporting Data also give other information including 
the camera tilt angle and azimuth, and the scale factor. The scale 
factor for the high-resolution photographs is based on the original size 
of the frame being approximately 55 mm across and each framelet about 
2.54 mm wide. It is understood that these measurements are fairly 
constant.* if the photograph at hand covers only part of a frame so 


* Verbal communication with Dr. Henry J. Moore. 



that measurements cannot be made from corner to corner, the above 
information makes it possible to determine the scale nevertheless. 

This is done simply by multiplying the given scale factor by the ratio 
of framelet width of photograph at hand to 2.54 mm. 

These methods were applied for determination of the scale of the 
photographs used in the analysis. 

Measurements of Boulder and Track 

The boulders and tracks were first measured on high- resolution 
Orbiter frames (approx. 39.8 cm wide) . Then the measurements were 
repeated for most of the boulders on photographs further enlarged, 
having a scale roughly five times that of the high-resolution frames. 
The two measurements were averaged and the resulting values of boulder 
diameter and track width were used in the analysis. Equal weight was 
given to both measurements because, although objects appeared larger on 
the enlarged photographs, their boundaries were more blurred and harder 
to define. The two measurements differed from each other by an average 
of about 10 per cent, and the maximum difference was about 30 per cent. 
Such a variation is not surprising, since the smaller boulders were 
close to the limit of resolution of the photography. J 

In the analysis of the boulder tracks it was assumed that the 
boulders were spherical. To make this assumption as valid as possible, 
only boulders appearing equidimensional on the photographs and leaving 
relatively smooth and well defined tracks were selected for study. 

Estimation of Slope Angle 

The slope angle was estimated from shadow relationships on the 
photograph, using the sun angle as listed in the supporting data. For 
example, if it can be assumed that a boulder is spherical, the slope 
can be calculated from measurements of boulder diameter, track width, 
and length of shadow cast by the boulder. This slope will, of course, 
be in the direction of the shadow which may not be the direction wanted. 
A better determination results from the shadow cast by a crater rim or 
a relatively horizontal upper surface adjacent to a rille. All these 


methods, however, depend on certain assumptions based on the investi- 
gator^ interpretation of the photograph and may therefore be consid- 
erably in error. These methods were used to estimate the slope angle, 
ct, for the boulders analyzed. Some slope angles were also provided by 
the USGS through their photogrammetry procedure. 

Material Properties Used 

A rock density of 2.7 gm/cc, a soil density of 1.6 gm/cc, and a 
cohesion of 1 x io 3 dynes/cm 2 were used in analyzing the selected 
boulder tracks.* 

RESULTS 

The results of the analysis are presented in Figures 2-9 and 2-10. 
Figure 2-9 shows the relationship between the track width to boulder 
diameter ratio and the friction angle. Figure 2-10 shows the frequency 
distribution of the results. As shown in the figure, the friction 
angles determined as described in this chapter were between 27 and 41 
degrees for most of the boulder tracks analyzed. The average friction 
angle was about 34 degrees. 

Some comment on the potential usefulness of the type of plot 
presented in Figure 2-9 is in order. A fairly well-defined single 
curve is shown, even though the data represent analyses for slope 
angles, (x, between 0 and 30 degrees. Therefore, such a curve appears 
suitable for a rough estimate of the friction angle of lunar soils for 
a variety of conditions from boulder-track measurements. Even without 
knowing the slope angle, the resulting friction angle would in the 
extreme case be possibly 5 degrees in error but usually only 1 to 3 
degrees from the correct value based on this theory. If a better theory 
could be developed, similar graphs could be prepared, and fairly 
reliable values of the actual friction angle of lunar soils could be 


* Measurements on samples returned by Apollo 11 and 12 (subsequent to 
these analyses) indicate a rock density of about 3.1 gm/cc. 








easily determined from the boulder-track relationship. 

Detailed data and results on each boulder track studied are 
presented in Chapter 7 (Table 7-1A, B, C) where the boulder tracks 
are reevaluated using an improved theory. 



Chapter 3. MODEL STUDIES OF THE FAILURE MECHANISM 
ASSOCIATED WITH A SPHERE ROLLING DOWN A SOIL SLOPE 


Model experiments were conducted to form a rational basis for 
theoretical developments, and to gain some understanding of the three- 
dimensional nature of the sphere-soil slope interaction phenomenon. This 
chapter describes how the models were constructed, presents photographs 
of the results, and discusses possible interpretation of the observed 
features. A failure mechanism is proposed, which may also be applicable 
to many cases of wheel-soil interaction. 

MATERIALS AND EQUIPMENT 

The models were constructed ih a rectangular pan approximately 
50 cm long, 30 cm wide, and 12 cm deep, using a mixture of ten parts fine 
to medium air-dry sand and two parts plaster of paris. Gradation of the 
sand (PCA Lapis Lustre #0 fine sand) and the mixture are shown in Figure 
3-1. Figure 3-1 also compares this gradation to the sand used for the 
rolling sphere tests at WES (Yuma Sand) , and to that of actual lunar soil 
(Apollo 11 sanples) . 

Additional test equipment consisted of: 

1. a glass cylinder or jar used for compaction, 

2. a shaker and a dropper used for deposition of intermittent 
dark layers and vertical markers, respectively, and 

3. a sphere (diameter of 12.10 cm, density of 3.69 gm/cc) . 

The sphere was rolled over the sand/plaster of paris models, which 
were then hardened by saturation. It was noted that the top of the models 
settled about 1 mm during saturation. After saturation and solidification, 
density measurements were made on pieces of the models. From these measure- 
ments, the following approximate values of density and void ratio were 
obtained: 


density before saturation = 1.59 gm/cc 
void ratio before saturation = 0.68 
density after saturation = 1.61 gm/cc 
void ratio after saturation = 0.66. 
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MODEL PREPARATION AND TESTING PROCEDURE 


The models were prepared and tested as follows: 

1. A 1-cm-thick layer of clean coarse sand was placed on the bottom 
of the pan to serve as a drain. 

2. The first layer (about 1-cm thick) of the sand and plaster of 
paris mixture was placed gently on the coarse sand layer and 
leveled off with the top of the pan. 

3. This layer was compacted by rolling a glass jar once over the 
surface. Additional material was then added and leveled to 
make the layer exactly 1-cm thick. 

4. A thin (approximately 1 mm) layer of a cement and coal powder 
mixture was deposited using a shaker. 

5. The next layer of sand/plaster of paris mixture was placed in 
the same fashion. The soil directly on the dark layer of cement 
and coal powder had to be placed with great care to prevent 
disturbance. The model was built up to the desired height 
using 1.0- and 0.5-cm-thick layers. 

6. Vertical markers were inserted to enable the determination of 
horizontal movements. This was accomplished by making a 1-mm- 
wide slot vertical to the drainage layer with a spatula. By 
moving the spatula gently back and forth it was possible to 
make the slot stay open when the spatula was removed. Coal 
powder was then deposited into the slot using a dropper. 

7. A template with small holes in a hexagonal pattern was placed 
on the surface of the model, and a cement-coal powder mixture 
was shaken over the template. Removal of the template left a 
hexagonal pattern of dark dots on the top surface of the model. 
Distortion of this pattern indicated surface movements. 

8. The sphere was placed on one end of the model, and the model 
was gradually tilted to an angle of about 13 degrees, at which 
point rolling was impending. The sphere was rolled the desired 
distance and taken off, and the model tilted back to a horizontal 
position. 
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9. The model was water-saturated slowly from the bottom? this 
process took about 1 hour. 

10. The model was left to harden and cure. 

11. After solidification, the model was sawed along desired sections. 
RESULTS 

Photographs of the models are presented in Figures 3-2 through 3-6. 
Deviations from straight lines in both top view and section photographs 
show the extent and nature of deformations. 

Figure 3-2 gives a top view of Model 1 and the numbering of the 
cross sections. Knife marks along the edges show where the model was 
subsequently cut to expose cross sections. One longitudinal section 
was made down the center of the track. Cross sections of Model 1 are 
shown on Figure 3-3, and the longitudinal section is shown on Figure 3-4. 

Model 2 was constructed primarily to study general shear in greater 
detail by using thinner layers (0.5 cm) and to study lateral movements. 

The cross sections of Model 2 are shown in Figure 3-5, a longitudinal 
section in Figure 3-6. 

DISCUSSION 

The most noticeable features on the preceeding photographs appear to 

be : 

1. Forward movement of soil as shown in both top-view and longitudi- 
nal section photographs. Lateral movement of soil is relatively 
small. 

2. Volume change of soil (compression directly under the sphere and 
dilation to the sides and front) as shown in the cross sections. 

3. Shear planes, particularly in the longitudinal sections, and 
also in the cross sections. 

Forward Movements 

Perhaps the most interesting observation from these models is the 
pronounced forward movement of soil immediately under the track. The nature 
of this movement is shown in Figures 3-2, 3-4, and 3-6. It is of interest 



surface deformation. The small 







3-3a. Cross sections of model 1 







3-3b. Cross sections of model 1 
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to consider the subsurface movements in more detail and to see if the 
subsurface and surface movements are related. 

Subsurface movements are best seen on the longitudinal section of 
Model 2 (Figure 3-6) where forward movements are represented by bending 
of the originally vertical dark markers. To illustrate certain details, 
the same section with some added lines is presented in Figure 3-7a. An 
enlargement of part of the model (Figure 3-7b) gives the position of 
the sphere when it stopped — just back of the fourth vertical marker. 

On the enlargement, the third vertical marker is superimposed on the 
location of the fourth. Two conclusions can be drawn: 

1. Most of the forward movement takes place before the edge of 
the sphere encounters the soil, as shown by the bending of 
the fourth marker. 

2. Some additional forward movement takes place as the point on 
the sphere where soil contact is made moves to the position 
vertically below the center of the sphere (indicated by the 
bending of the third marker) . 

The net effect of rolling a sphere over a soil slope is, therefore, 
a forward bending of the originally vertical sections. It seems 
that the bent shape could be reproduced simply by bending the section 
forward as if it were a steel spring. Minor stretching would be 
required, and it is believed that this is associated with general shear 
close to the surface. It should be noted that, while bending of these 
sections (markers) extends to a depth of approximately 5 cm, no general 
shear planes can be seen below a depth of 1.5 cm from the bottom of the 
track. Hence, subsurface deformations observed must be largely due to 
shearing distortion and volume change. 

To compare forward movements as measured from the longitudinal section 
of Model 2 with those measured from the top-view photograph of Model 1, it 
is desirable to identify the movements on the top-view photograph of Model 1 
more clearly. Figure 3-8 is the same top-view of Model 1 as that in Figure 
3-2, with grid lines added. The line around the track outlines the approxi- 
mate outer limit of surface deformation, and the inner closed line is the 
approximate limit of soil-sphere contact. In this case, a track width that 
is slightly smaller than the crest- to-crest distance is defined. 
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A comparison of forward movements as measured on Figure 3-8 and from 
the longitudinal section of Model 2 (Figure 3-7a) is presented on Figure 
3-9. This figure shows that the movements are comparable and that, for 
this material, surface deformations are indicative of subsurface movements. 
Some deviation is seen at small track width or at the beginning of the 
track. This deviation is probably due to differences in the formation of 
the tracks for the two models. 

Volume Changes 

By studying the cross section photographs, it appears possible to 
distinguish zones of compression, dilation, and shear. In the following 
discussion these zones are determined, and the cross sections are analyzed 
in greater detail. 

Zones of compression were determined by noting at what points the 
distance between the dark interfaces decreased from the original values. 
(Linear measurements were made using a magnifying glass with a built-in 
scale reading to one-tenth of a millimeter.) Zones of dilation were 
determined similarily by noting at what point the distance between two 
dark interfaces increased from the original values. These criteria are 
illustrated in the cross sections in Figures 3-10 and 3-11. The region 
directly below the sphere is a zone of volume decrease. Regions immediately 
under the crests of the track are zones of dilation. The region in front 
of the sphere is also a zone of dilation (see Figure 3-11, cross section 4, 
and photographs of the longitudinal sections) . 

The magnitude of volume changes was measured from the photographs 
using a planimeter. Compression or volume decrease is represented by the 
segment of the circle below the respective reference line drawn parallel 
to the original surface or dark interface (Figures 3-10 and 3-11) . Dila- 
tion or volume increase is represented by the area above. the reference 
line and is associated with the crests of the track. The results of these 
measurements, as given in Figure 3-12, present dilation and compression 
versus the depth below the original surface of Model 1. These results 
are compared below with those of pressure under a sphere as a function of 
depth . 
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An assumption in many problems is that deformations are proportional 
to the applied pressure. This assumption is generally considered valid 
for sands, although the proportionality coefficient, k g (Terzaghi, 1943), 
is a function of depth. It is of interest to determine to what extent 
this proportionality is valid for the case of the rolling sphere. This 
can be accomplished by considering Boussinesq ' s solution of the pressure 
distribution beneath a circular loaded area. The pressure distribution 
at the soil-sphere contact is believed to be approximately parabolic. 

In the vicinity of the second cross section of Model 1 the track width 
was approximately 8 cm. Since the actual contact area is somewhere 
between circular and semicircular (Hovland and Mitchell, 1970), an average 
width of 6 cm was assumed for comparison purposes. Figure 3-13 compares 
pressure and deformation as functions of depth. The deformation curves 
are actual measurements from the cross sections. It is apparent that the 
pressure curve and the deformation curve for the second cross section are 
almost identical. These curves should also be compared with Figure 3-12. 
From the similarity thus established, it can be concluded that the de- 
formations are approximately proportional to the applied pressure. 

General Shear 

General shear has been shown to develop predominately in the forward 
direction (Model 2, Figure 3-la ) . The lateral shear in the cross section 
of Model 2 is probably a consequence of forward shearing. 

To explain some of the details observed on the longitudinal sections 
associated with general shear, the sequence of rolling and shearing of a 
sphere is diagrammed in Figure 3-14; (a) the soil before the sphere has 
rolled over it; (b) the shear planes and deformations that develop as 
the sphere rolls to the right, to position B; and (c) the appearance of 
the section after the sphere has rolled past. These same -features can be 
observed in the photographs. 

Observations of particular interest are: 

1. Shear planes do not appear to originate at the sphere surface 
(see Figure 3-1 4b) , as is often considered the case with shear 
under a loaded footing. ' 
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2. At the time a shear plane develops, its lower end is roughly 
parallel to the sphere surface. It diverges from a direction 
parallel to the sphere surface as it proceeds up and forward. 

3. A wedge of soil is apparently pushed up and forward. 

4. Movement along any one shear plane is restricted, as implied 
by the relatively short but constant displacements. 

5. Shear planes are spaced at relatively constant intervals as 
can be seen from the spacing of the teeth marks (compressed 
remnants of once active shear planes) . 

It is possible to conclude that the shearing process is not continuous, 
but consists of separate small shear failures spaced at relatively constant 
intervals with respect to space and time. The spacing of the shear planes 
as well as the amount of movement along any one shear plane is undoubtedly 
a function of soil type. Hence, although the rolling of a sphere down a 
slope may appear continuous to the naked eye, it is possible that the 
acceleration of the sphere changes slightly from one shear plane to the 
next. 

The fact that shearing along any one shear plane appears to be 
restricted to a relatively small amount is perhaps associated with the 
change in direction of the shear plane and of the stress causing the 
movement. Most of the movement along a shear plane probably takes place 
when the plane is first formed. At that time, it is directed forward at 
the critical angle of 45° — <f>/2 to the direction of the major principal 
stress (Scott, 1963) . Hence, one can also estimate the direction of the 
principal stresses along the shear plane, provided it is in the position 
where it was first formed. As the sphere moves forward, the shear plane 
bends or turns toward a more vertical position. As the sphere passes, 
the shear plane is again bent down toward a final, more horizontal position 
(Figure 3-14) . 
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PROPOSED FAILURE MECHANISM 

f The proposed hypothesis of a failure mechanism presented below is 

restricted to a rigid sphere rolling on a deformable material such as 
a sand slope. For such a material, the failure mechanism is dependent 
on the state of compaction. In a very loose state, the material would 
’ be compressible and the deformations would consist of both volume change 

and shearing distortion, but no general shear planes would develop. In 
a very dense state, general shear would be more significant, although 
volume change and shearing distortion would also take place. The usual 
case would involve all three phenomena: Initial volume change accompanied 

by shearing distortion would be followed by general shear if the material 
i were sufficiently dense. 

The following steps appear to be involved as a sphere rolls over 
a sand slope: 

1. Compression occurs under the sphere. 

2. Initially vertical sections are bent forward. 

3 . A horseshoe-shaped zone of dilation , which extends from the 
sides and around the front of the sphere, develops. This zone 
moves forward with the sphere (analogous to a congressional 
wave) . The external evidence of this zone of dilation is the 
bow wave that forms in front of the sphere as material accumu- 
lates . 

4. When volume changes and shearing distortion can no longer account 
for all the soil displacement, general shear planes develop. The 
development of general shear planes is probably influenced by the 
magnitude of shearing distortion and the magnitude and direction 
of the major principal stress. The amount of shearing distortion 
that can occur without the development of shear planes is propor- 

, tional to the confining pressure. Confining pressure increases 

with depth, and the applied stress at the sphere-soil contact 
increases approximately parabolically from the edge of sphere-soil 
contact. Therefore, shear planes would develop where the deforma- 
tion reaches a critical magnitude with respect to the confining 
pressure. These planes would be oriented approximately in a 
direction of 45° — 4>/2 to the direction of the major principal 
stress. 
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5. Sliding along the shear planes continues only for a short time 
and distance, as previously described. 

6. After a certain distance, a new shear plane develops and the 
shearing cycle repeats itself. 

7. As the sphere rolls down the slope, deformations assume a final 
position, as shown on the photographs, and excess material from 
the bow wave is wasted to the sides to form the crests of the 
track. 

CONCLUSIONS AND RECOMMENDATIONS 

Models such as those used for this study have proved valuable in the 
investigation of surface and subsurface deformations. The vertical dark 
markers proved to be particularly helpful. 

Due to the limited model length, the track was short. It is possible 
that some additional effects would develop with long, continuous rolling. 
The proposed failure mechanism is reasonable . 

Resistance to rolling on sand appears to be caused by volume change, 
shearing distortion, and general shear. At higher velocities, inertia 
effects may also be important. Determination of the relative contributions 
of these phenomena to the total resistance requires further investigation. 
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Chapter 4. PRESSURE DISTRIBUTION AND SOIL FAILURE 
BENEATH A SPHERICAL WHEEL IN AIR-DRY SAND 

To develop a theoretical explanation of the rolling sphere-soil slope 
interaction phenomenon, the pressure distribution beneath the sphere, in 
addition to the failure mechanism, must be understood. A failure mech- 
anism associated with a sphere rolling down a soil slope was proposed in 
Chapter 3. This chapter describes studies to determine the pressure 
distribution. 

Experiments were carried out using an instrumented spherical wheel 
at the U. S. Army Engineer Waterways Experiment Station (WES) , Vicksburg, 
Mississippi. The assistance provided by the WES to this effort was funded 
through NASA Defense Purchase Request (DPR) H-58504A from MSFC to WES. 

The results provided information on the following: 

1. The distribution of contact pressure beneath a towed spherical 
wheel rolling on air-dry sand, arid 

2. The appearance of the resulting track, revealing surface 
features of soil deformation and failure. 

The equipment, facilities at WES, and soil used are described first. 
This is followed by a description of the testing procedure and results. 

The discussion and analyses are primarily concerned with the contact 
pressure distribution. 

EQUIPMENT, FACILITIES, AND SOIL 

Equipment 

A special spherical wheel was designed and instrumented in order to 
measure normal pressure at the sphere-soil interface during rolling. 

The spherical wheel was intended to simulate a solid spherical boulder. 
Therefore, it was designed to act as a rigid body. Figures 4-1 and 4-2 
show side and front views of the spherical wheel. The outer shell was 
1/2 inch thick. Deformations were calculated using an axisymmetric 
finite element program based on shell theory (M. Khojasteh-Bakht , 1967). 
For a 1/2-inch- thick hemispherical shell, 24 inches in diametef , loaded 
by an axisymmetric parabolic pressure distribution, the maximum deflection 
of the shell, which takes place directly under the point of maximum 

found to be 0.0011 inches for a maximum pressure of 40 psi. 


pressure, was 
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MATERIAL ’• ALUMINUM 


Fig. 4-1. Side view of spherical wheel. 
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A row of semiconductor diaphragm gages was located along the 
periphery of the spherical wheel as shown in Figures 4-1 and 4-2. The 
gages were designed to measure the normal pressure at the sphere-soil 
contact as the sphere rolled through the soil. A total of seven gages 
was used; their location and numbering are shown in Figure 4-2. Addi- 
tional positions (holes) were provided to make it possible to relocate 
the gages for better coverage of pressure at large sinkage. Such 
relocation turned out to be unnecessary, and the holes were filled using 
a screw. 

Figure 4—3 shows details of a diaphragm gage. The gages were 
screwed in from the inside with a special key and locked in position 
with a lock nut as shown in Figure 4-3. The diaphragm was designed 
using formulas presented by Timoshenko and Goodier (1951) . Theory 
indicated that for a 40-psi normal pressure applied to a 1/2-inch-diameter 
and 0.012-inch- thick diaphragm of cold rolled 1045 steel, the maximum 
strain would be about 700 microinches per inch for a maximum deflection 
of about 500 microinches. This design proved to be satisfactory. 

The semiconductor strain gages used are described by the manufac- 
turer (BLH Electronics, Inc.) as follows: 

Type: SPB 3-12-12 

Resistance: 120 ohms 

Gage Factor: +116 

Length: 0.12 inches . 

The gages were calibrated at WES. Wires from the semiconductor 
gages were run to a slip-ring assembly, and resistors were added to 
complete a Wheatstone bridge. The slip-ring assembly, a panel of added 
resistors, and the spherical wheel in position in the carriage are 
shown in Figure 4-4. To prevent diaphragm damage by point loads from larger 
soil grains, a strip of thin rubber was stretched to cover each gage 
face. Lastly, fine uniform sand was glued to the spherical wheel 
surface to make the surface rough like the surface of a boulder. 

Facilities 

The spherical wheel was towed on the soil car shown in Figure 4-5. 

This car was 64 inches wide, 36 inches deep, and more than 50 feet long. 




4-4. Spherical wheel in carriage 
(photograph also shows the 
slip-ring assembly and panel 
of added resistors). 




Camera 


Drawbar 


Spherical wheel 


Fig. 4-5. Testing facilities at WES 



Figure 4-5 also shows the cable which applied the towing force to the 
carriage (upper right) , and the movie camera (left of the floodlight) 
which followed the carriage and filmed the sphere rolling through the 
soil. 

The pressure registered by the gages, as well as the towing force, 
wheel load, velocity, acceleration, and inertia effects of the carriage 
were recorded continuously on an oscillograph as the spherical wheel 
rolled through the soil. 

Soil 

A gradation curve for the Yuma sand used in these tests is shown in 
Figure 3-1. Yuma sand is a uniform very fine sand with an effective 
grain size of about 0.08 mm and a uniformity coefficient of about 1.5. 

In its loosest state (D^ = 0%) , air-dry Yuma sand has a friction angle 
of about 35 degrees. In its densest state (D^ = 100%), air— dry Yuma sand 
has a friction angle of about 47 degrees. The air-dry soil densities 
corresponding to the loosest and densest states are 1.39 gm/cc and 
1.66 gm/cc, respectively. The specific gravity of Yuma sand is 2.67. 

TESTING PROCEDURE 

Testing proceeded through the following steps: 

1. The soil in the car was prepared to the desired density using 
techniques established at WES as described in Appendix A. 

With this procedure, a particular density could be repeated to 
within plus or minus 2 per cent. The Yuma sand used in this 
study was compacted to give a cone index gradient value of 
approximately 2.00 meganewtons per cubic meter (7.42 lb/in 3 ). 

2. The pressure gages were calibrated, and the wheel load to be 
used was selected. 

3. The spherical wheel was rolled approximately two revolutions 
on the soil car at a velocity of about 2 feet per second. 

4. The pressures registered by the gages and other parameters of 
carriage motion were recorded. 

5. The tracks formed by the rolling wheel were measured and 
photographed. 
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RESULTS 

The spherical wheel was tested under five wheel loads. These wheel 

loads , although applied at the axle r gave an effect identical to solid 

spheres of different densities r since the sphere was rigid. The loads 

were selected to give equivalent density ratios (density of rock/density 

of soil = Y /y ) corresponding to solid spheres built to investigate the 
r s 

dynamics of spheres rolling down a sand slope (see Chapter 6) . The loads 
used and the corresponding equivalent density ratios are shown in Table 4-1. 


Table 4-1. Wheel loads. 


Wheel load 
N 

Corresponding 
density, Y r 

gm/cc 

Soil density 
Y s 

gm/cc 

Density ratio 

y h 

'r s 

635 

0.55 

1.57 

0.35 

1241 

1.08 

1.58 

0.68 

2100 

1.83 

1.58 

1.16 

2503 • 

V.. • 

2.18 

1.57 

1.39 

4275 s : 

3.71 

1.57 

2.36 


The resulting tracks showing soil deformation and failure under these 
wheel loads, and the measured distribution of contact pressure as pro- 
jected on a horizontal plane are shown in Figures 4-6 through 4-15. 

The pressures shown in Figures 4-6 through 4-15 were plotted from the 
oscillograph data presented in Appendix C. 

DISCUSSION 

In general, results show that the contact pressure normal to the 
surface of the spherical wheel has a bell-shaped distribution. In the 
discussion below, the pressures, the resulting tracks, and some analytical 
formulations are considered. 


a) Track resulting from a wheel load of 635 newtons. 


Fig. 4- 


vert. scale: 1" = 10 cm 
hor. scale: 1" = 10 in. 


cross section, St. 29.0 



longitudinal section, St. 29.55 
b) Sections for the above track. 

6. Track and sections for a wheel load of 635 newtons (143 lb) 



SPHERE MOVED TO RIO NT 
PRESSURE IN KN/lR* 



SCAl E : 


PRESSURE 





a) Track resulting from a wheel load of 1241 newtons. 


vert, scale: 1" = 10 cm 
hor. scale: 1" = 10 in. 


cross section, St. 25.0 



longitudinal section, St. 25.0 
b) Sections for the above track. 


Fig. 4-8. Track and sections for a wheel load of 1241 newtons (280 lb). 
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I 

I a) Track resulting from a wheel load of 2100 newtons. 


vert, scale: 1" = 10 cm 
hor. scale: 1" = 10 in. 



cross section, St. 19.5 


1 _ 

z = sinkage 

I 

longitudinal section, St. 19.5 
b) Sections for the above track. 



Fig. 4-10. Track and sections for a wheel load of 2100 newtons (472 lb). 






a) Track resulting from a wheel load of 2503 newtons 

vert, scale: 1 " = 10 cm 
hor. scale: 1" = 10 in. 


cross section, St. 19.0 


longitudinal section, St. 20.7 
b) Sections for the above track. 


Fig. 4-12. Track and sections for a wheel load of 2503 newtons (564 lb) 








a) Track resulting from a wheel load of 4275 newtons 


vert, scale: 
hor. scale: 


1" = 10 cm 
1" = 10 in. 


V" original contact 
^ surface 



caved 


\\ surface 



cross section, St. 24.1 


V ^ caved surface 
(angle of repose) 









Pressure at the Sphere-Soil Interface 

The contact pressure normal to the surface for each wheel load was 
projected or plotted on a horizontal plane. Figure 4-16 shows the 
contact pressure plotted on the curved surface of the sphere along a 
longitudinal section through the center of the track. Figure 4-16a 
shows the actual location of the pressure distribution for the lightest 
wheel load used (635 newtons) , and Figure 4-16b shows the same infor- 
mation for the heaviest wheel load (4275 newtons) . The figure also shows 
the direction of the resultant of towing force and wheel load, and the 
direction of maximum pressure. 

The ratio of towing force to wheel load increased with increasing 
wheel load, as shown by Figure 4-17. (For these wheel loads, slip varied 
as shown in Figure 4-20.) 

As illustrated in Figure 4-16, the direction of the resultant of 

towing force and wheel load is not quite parallel to the direction of 

maximum pressure. This deviation was investigated and is plotted on 

Figure 4-18. (3 is the angle between the vertical and the direction 

s 

of maximum pressure.) Since the resultant of towing force and wheel 
load must be equal and opposite to the resultant of the sum of pressure, 
the resultant of the sum of pressure must act slightly below or behind 
the location of maximum pressure. This conclusion is reasonable since 
the contact area is considerably wider toward the back. 

To check the reliability of the measurements, the sum of all 
pressure within the pressure distribution was computed by determining 
areas between the pressure contours in Figures 4-7, 4-9, 4-11, 4-13, 
and 4-15, using a planimeter, and multiplying each area by its average 
pressure. To add these incremental forces (pressure times area) , however, 
their components parallel to the direction of the maximum pressure had 
to be found. This was done by dividing each force by the cosine of the 
angle between the direction of each respective force and the maximum 
pressure. Table 4-2 shows a comparison between the resultant of wheel 
load and towing force and the sum of pressure (reaction force) computed 
as described above. 
















Table 4-2. Comparison between applied and measured loads 


Wheel load 
(N) 

Towing force 
(N) 

i 

Resultant 

(N) 

Reaction force 
computed from 
pressure diagrams 
(N) 

635 

79 

640 

611 

1241 

205 

1255 

1216 

2100 

520 

2165 

2178 

2503 

723 

2610 

2546 

4275 

1850 ' 

4650 

4637 


There is good agreement between the resultant and the reaction force as 
presented in Table 4-2. The average deviation is about 1 per cent, and 
the maximum deviation is about 4 per cent. On this basis, we can 
conclude that the pressure gages reliably measured the normal contact 
pressure between the spherical wheel and the soil. 

While the sum of pressure (reaction force) within the bell-shaped 
distribution increased in proportion to the applied loads, the maximum 
pressure did not. The maximum pressure resulting from the increasing 
applied loads is plotted in Figure 4-19. This plot implies that the 
maximum pressure increases until the soil begins to fail. For higher 
loads, the maximum pressure remains nearly constant, and added resis- 
tance is supplied by greater sinkage of the sphere, resulting in a 
larger contact area. It appears reasonable to conclude that shear 
failure is more directly a function of the maximum pressure than of the 
average bearing pressure. 

Slip (defined as the distance revolved minus the distance traveled 
divided by the distance traveled) was found to increase with increasing 
wheel loads or density ratios, as shown in Figure 4-20. All slip was 
found to be negative. The negative slip does not necessarily mean that 
there is a net shear stress acting on the sphere surface. As shown from 
the model studies (Chapter 3), a soil wedge is pushed up and forward, 
which would cause an opposite shear stress on the front part of the 
sphere. This interpretation agrees with the previous experimentation 
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RESULTANT OF WHEEL LOAD AND TOWING FORCE IN NEWTONS 


Fig. 4-19. Maximum pressure vs resultant of Wheel load 

and towing force • 



Fig. 4-20. Slip vs density ratio 



with towed wheels (Green and Murphy, 1965) . 


Tracks 

Photographs of the tracks (Figures 4-6 through 4-14) show certain 
repeated features that agree with the failure mechanism proposed in 
Chapter 3. Other visible features, important to an understanding of 
wheel-soil interaction ares 

1. The track for the lightest wheel load (Figure 4-6a) shows no 
evidence of general shear, which agrees with the above inter- 
pretation of Figure 4-19. 

2. The tracks for the lighter wheel loads (Figures 4-6a and 4-8a) 
also show ripples in the center of the track. Such ripples are 
not seen in the tracks resulting from the heavier wheel loads. 

3. Two sets of well-defined shear planes are clearly shown in the 
tracks resulting from the heavier wheel loads (Figures 4- 10a, 
4-12a, 4-14a) . One set runs in front of the spherical wheel 
across the track with a slight curvature. The other runs 
approximately perpendicular to the first set, diverging out 
from the direction of travel at an angle of about 30 degrees. 

4. Very little lateral movement of soil took place, as shown by 
the fact that the small longitudinal grooves seen in Figures 4-6 
through 4-15 were not distorted. (These grooves were left on 
the soil surface by the screeding or leveling operation.) This 
agrees with the results of model studies described in Chapter 3. 

Experimental information on shear zones (the zone within which 
failure planes could be detected in the soil) and the number of failure 
planes within the shear zone are presented in Figure 4-21. The maximum 
width of the shear zones extending both forward and to the sides was 
measured from the crest of the track for each track; the results are 
plotted (Figure 4-21a) as a function of the density ratio. This figure 
should provide a valuable comparison to theoretically predicted shear 
zones (example. Chapter 5, p. 94). The number of failure planes within 
the shear zones were counted and are presented (Figure 4-21b) as functions 
of track depth. This plot suggests that there were between one and two 
failure planes per centimeter depth of soil within the shear zones. Most 




RATIO OF MAXIMUM WIDTH OF SHEAR 
ZONE TO AVERAGE TRACK WIDTH 


a) Density ratio vs width of shear zone. 



b) Number of failure planes within the forward shear zone 
as a function of track depth. 


Fig. 4-21. Experimental information on shear zones 
and failure planes within the shear zone. 
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of the data for Figure 4-21 are from the spherical wheel rolling on Yuma 
sand. It is interesting that the one point from the models reported in 
Chapter 3 fits quite well within the data, as shown in Figure 4-21, even 
though the models were constructed of a somewhat different material. 

Analytical Considerations 

It is possible to approximate the bell-shaped pressure distribution 
with a parabolic function, as shown in Figure 4-22. Therefore, we may 
represent the contact pressure in a direction parallel to the direction 
of the resultant of wheel load and towing force by the expression 

.2 


0 = 0 - ay' 

max 1 


(4-1) 

is the radius of 


When y = y , 0=0. Therefore, a = O /y 2 ; v 

mx max max' 

a circular area equal to the actual approximately semicircular contact 

area* Therefore, 


Try 


max 


2L 

2 4 ' 


(4-2) 


where w equals the crest-to-crest track width. From this egression. 


1 2 
W 


Then 


8 


max 


w 


(4-3) 


(4-4) 


The total applied force is the resultant of wheel load and towing 
force. The reaction to this resultant is the sum of all pressure within 
the pressure distribution described by Equation (4-1) , which acts in a 
direction opposite to the total applied force. (It is assumed that net 
shear along the sphere surface is zero.) The total applied force and 
its reaction can be equated by the solid of revolution approach to give 


* - * f “ f 2 »>« -l [""* ( CT max - o)*> . (4-5. 

J o J o 

V 

where R is the resultant of wheel load and towing force. This approach 
assumes that the actual approximately semicircular area can be represented 
by an equal circular area. Substituting the value for "a" from Equation 
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l 


(4-4) into Equation (4-5), and carrying out the integration gives 


R 



a 

max 


(4-6) 


or since the area equals A 


TT w 
2 4 


R = ! 


a a 


max 


(4-7) 


from which 


a 

max 



(4-8) 


where q * R/A is the average pressure often referred to as "unit bearing 
pressure." Equation (4-8) shows that the maximum pressure beneath the 
spherical wheel is about twice the average unit bearing pressure. 

To determine how closely the above equations describe the contact 

pressure beneath the spherical wheel, the resultant, R, was computed 

using Equation (4-6) with measured a , and the maximum pressure was 

max 

computed using Equation (4-8) with measured R. Table 4-3 presents 
the comparison. 


Table 4-3. Comparison between applied and 
calculated loads and pressures. 


Wheel load 
N 

Applied R 
N 

Calculated R 
N 

Measured a 

max 

KN/m 2 

Calculated a 

max 

KN/m 2 

635 

640 

630 

73 

74 

1241 

1255 

1330 

93 

88 

2100 

2165 

2170 

102 

102 

2503 

2610 

2680 

112 

109 

4275 

4650 

4790 

106 

103 


The comparison shows a very good agreement, and it is believed that 
for a spherical wheel in air-dry sand, the ratio between a 

max 

about 2. 


and q is 
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Combining Equations (4-1) , (4-4) , and (4-8) leads to an equation 
for the contact pressure beneath a sphere in air-dry sand, 

. 2 - 


O - 2q 


MiO]' 


(4-9) 


CONCLUSIONS 

It has been illustrated in this chapter that the distribution of 
contact pressure beneath a rigid sphere or spherical wheel in air-dry sand 
is bell-shaped. This bell-shaped distribution can be closely approxi- 
mated by an equivalent parabolic function. 

Shear failure appears to be more directly related to the maximum 
pressure beneath the sphere than to the average pressure or unit bearing 
pressure. The maximum pressure beneath a rigid sphere or spherical 
wheel in air-dry sand is shown to be about two times the average pressure. 

The pattern of deformation and failure of the soil appears to agree 
with the failure mechanism proposed in Chapter 3. 

It is believed that the results clarify the behavior of a rigid 
sphere rolling in air-dry sand. Further theoretical analysis can be 
based upon these findings. 
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Chapter 5. THEORETICAL STUDIES 

Investigation of the rolling sphere-soil slope interaction problem 
has been emphasized only recently. To the writer's knowledge, the only 
existing theory is that proposed by Lin and Haythornthwaite (1969) . They 
studied ping-pong-ball size spheres rolling down a slope of granular media 
at constant velocity. They proposed a solution based on the assumption 
that failure of the soil in front of the sphere is caused by bulldozing. 
Their theory and experimental data are for a cohesionless (c = 0) soil. 

Existing theories for wheel-soil interaction have been either quasi- 
theoretical (Bekker, 1956), or empirical (Reece, 1965-1966; Roth, 1960; 
Waterways Experiment Station, 1954). Due to mathematical complexities 
associated with three-dimensional analyses, it has also been suggested 
that the wheel-soil interaction problem be assumed to be two-dimensional 
(Yong and Osier, 1966) . Model studies of soil deformation under wheels 
have been carried out to investigate two-dimensional behavior (Wilson and 
Krzywicki, 1966) . 

Due to the quasi-theoretical or empirical nature of existing theories, 
much design has depended on experimentation for each new specific condition. 
It has not been demonstrated to what extent the failure mechanism observed 
in a two-dimensional model approximates the three-dimensional wheel-soil 
or sphere-soil interaction. 

This chapter begins with a summary of the significant conclusions 
from the results presented previously as a basis for theoretical develop- 
ments. This is followed by step-by-step development of a theory to 
provide a method by which soil properties can be evaluated from boulder 
track data. Simplifying assumptions are used in order to make the 
problem tractable. 

BASIS FOR THEORETICAL DEVELOPMENT 

From the results of the investigations of the failure mechanism 
using model studies and the pressure distribution of the soil-sphere 
contact using an instrumented spherical wheel, the following 



characteristics of behavior have been determined: 

1. Resistance to rolling appears to be a function of volume change, 
shearing distortion, general shear, and soil inertia. 

2. For loose sands at the beginning of rolling, soil deformations 
consist of volume change (compression) and shearing distortion. 

3. When volume changes and shearing distortion can no longer account 
for all the soil that must be displaced, general shear planes 
develop. These shear planes probably develop where the deforma- 
tions reach a critical magnitude with respect to the confining 
pressure. These shear planes are initially oriented approxi- 
mately in a direction of 45° - (f>/2 to the direction of the 
major principal stress. 

4 . It appears that a wedge of soil , originating approximately 
radially below the center of the sphere-soil contact area, is 
pushed up and forward during the shearing process. 

5. Initially vertical sections are bent forward as the sphere rolls 
down the slope and deformations assume a final position as shown 
in Chapter 3 . Material from the bow wave is wasted to the 
sides to form the crests of the track. 

6. Deformations (compression) beneath the sphere at any depth are 
almost directly proportional to the pressure at that depth. 

7. The distribution of contact pressure beneath a rigid sphere in 
air-dry sand is bell-shaped. This bell-shaped distribution 

can be closely approximated by an equivalent parabolic function, 

0 = 2q|l - (5-1) 

where 

a = contact pressure at any point parallel to the 
resultant of applied forces , R. 

q = average contact pressure parallel to 0, or 
R/contact area. 



y - distance from location of maximum pressure 
perpendicular to the direction of maximum 
pressure . 

w = crest-to-crest tract width. 
y max is related to w according to y^^ = w/2.825. 

8. General shear failure appears to be more directly related to 

the maximum pressure beneath the sphere than to the average 

pressure. The maximum pressure, a , is about 2q or* twice 

max 

the average pressure. 

9. The maximum width of the shear zone measured from the edge or 
crest of the track was found to be about 1.3 w in the forward 
direction for the heaviest sphere load and a sphere-to-soil 
density ratio of 2.2 for Yuma sand. The ratio of side shear 
to forward shear was found to be about 0.38 for all wheel 
loads . 

10. Negative slip increases with wheel load or sphere density, or 
with increased sinkage. 

Any theory must consider the above stated characteristics of 
behavior. Further requirements of theory are: 

1. Dynamic equilibrium must be satisfied. 

2. The three-dimensional nature of the problem must be recognized 
and satisfied or closely approximated. 

3. Soil and sphere properties must be accounted for. 

4. Soil deformation and volume change must be considered. 

5. General shear must be evaluated. 


6 . 


Inertia effects must be evaluated. 
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DYNAMIC EQUILIBRIUM 

A sphere rolling down a deformable surface, such as a soil slope, 
is shown in Figure 5-1. 


Fig. 5-1. Sphere rolling down a soil slope. 


The symbols in Figure 5-1 are defined as follows: 

W = weight of sphere. 

m *» mass of sphere = W/g. 

F = resultant soil reaction force. 

F^ = component of F parallel to slope. 

F^ = component of F normal to slope. 

r = radius of sphere. 

r = distance from center of sphere to 
F 

line of action of F. 

A 

I * moment of inertia of sphere =(2/5)mr . 



a = linear acceleration of sphere . 
u = angular acceleration of sphere . 
x,y = coordinate directions . 
a = slope angle . 

Negative slip will usually take place as a sphere rolls down a slope. 

Therefore, with respect to translation down the slope, a f ur. Any 

shear stresses that act on the sphere surface are assumed to be included 

as a part of forces F, F , and F . 

x y 

In applying Newton's Second Law using the method of dynamic equilib- 
rium, ma and Iu are considered acting opposite to their actual sense, as 
shown by the dotted arrows in Figure 5-1. We then have: 



= 0 = F - W cos a 

y 

- 0 - W sin a - F 

x 

= 0 = Fr„ - Iu, 

F ' 


ma. 


F = W cos a 

y 


F = W sin a - W - 
x g 


„ 2 „ r u 

F = — W — — . 

5 r „ g 
F 


(5-2) 


(5-3) 

(5-4) 


Equation (5-2) states that the normal component of soil reaction is 
always equal to the component of the weight of the sphere in the same 
direction. One cannot, however, conclude that the normal component, 

F^ , is independent of sphere motion, since it has been observed that at 
high velocities the sphere may even start bouncing. Equation (5-3) 
states that if the sphere is to be accelerating, the component of soil 
reaction parallel to the slope, , must be less than the component of 
the weight of the sphere in the same direction. In other words, if 
a = 0 (sphere resting on the slope or rolling with constant velocity) , 
the resultant of soil reaction, F, is exactly equal to the weight of 
the sphere, W, and has the same line of action. If the sphere is 
accelerating, the inclination of F moves toward the normal to the slope. 

Interpretation of Equation (5-3) is aided by Figure 5-2, which 
shows that: 

1. Positive values of F^/W represent soil resistance parallel to 
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the slope and opposite to the direction of motion of a freely 

rolling sphere. 

Here: 

a. The condition a = 90°, F x / W = If represents a sphere 
suspended in the air or falling at constant velocity. 

b. The condition a = 90° , F^/W - 0, represents a freely 
falling body. 

c. The positive vertical axis, F ^/W , represents necessary 
soil resistance for constant velocity rolling for a 
sphere in motion, or the force necessary to prevent 
rolling for a sphere at rest. 

2. Negative values of F^/W represent pull that must be applied 

to cause further acceleration. 

Here: 

a. The condition a = 0, F x / W ~ 0/ represents a sphere 
resting or moving at constant velocity on a horizontal 
surface. 

b. Any pull applied above and beyond the soil resistance 
will cause the sphere to accelerate. 

The resultant, F, in terms of F^ and F^, is 

F 2 = F 2 + F 2 . (5-5) 

x y 


Substituting F and F^ from Equations (5-2) and (5-3) into Equation 
(5-5) gives 



(5-6) 


The remaining unknowns in Equations (5-4) and (5-6) are F, a, u, 
and r . From experiments performed at WES, a and u are determinable 
for most spheres rolled. Therefore, F and r^ could be determined. 
This would give the magnitude and direction of the resultant soil 
reaction, F, as well as its point of action on the sphere surface. 
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While this relationship is useful for checking theory and experimental 
work, it could not be used for the evaluation of lunar boulder tracks 
since accelerations a and u are unknown. 

The amount by which the resultant, F, differs from the weight of 
the sphere, W, and the amount by which F deviates from the vertical is 
important to know. Figure 5-3 shows a plot of Equation (5-6) giving the 
value of F/W as a function of a/g. It may be noted that for a slope 
angle of 15 degrees, which was the average for 69 lunar boulder tracks 
studied (Hovland and Mitchell, 1970), F/W is between 0.97 and 1.00 
for reasonable values of a/g (a/g = 0 to 0.5). For other possible slope 
angles, F/W will vary between 0.85 and 1.10. Therefore, little error 
should be associated with the assumption that the resultant equals the 
weight of the sphere. 

Figure 5-4 shows the deviation of F from the vertical as a function 

of a/g for various slope angles. This figure shows that for a slope 

angle of 15 degrees and a/g approximately equal to 0.26, F will deviate 

from the vertical by about 15 degrees. This means that F is normal to 

the slope, and that F = F , and F =0. This situation (F =0) would 

y x x 

only arise for a frictionless slope. For a soil slope, experiments (to 
be reported in Chapter 6) indicate that the maximum value of a/g will 
be approximately 0.5 times the values for a frictionless slope. These 
maximum values of a/g were measured for light spheres rolling on dense 
moist soil. All other sphere and soil combinations would result in 
smaller values of a/g. The range of a/g for likely sphere-soil combi- 
nations would be from zero to 0.15. Therefore, the direction of the 
resultant would probably deviate only between 0 and 8 degrees from the 
vertical . 

In summary, the following conclusions can be added to the previous 
list of behavior characteristics: 

1. The magnitude of the resultant of soil reaction will usually 
be nearly equal to the weight of the sphere. 

. The direction of the resultant of soil reaction will usually 
be nearly vertical. 


2 
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The total soil reaction consists of resistance caused by volume 
change and shearing distortion, general shear, and inertia effects. This 
may be expressed per unit area by 

(5-7) 

where is resistance caused by volume change and shearing distortion, 

*lg resistance caused by general shear , and q^ is resistance caused by 
soil inertia effects. The contribution of these terms will be investi- 
gated in the following sections. 

In the above discussion and in the sections to follow, volume 
changes and shearing distortions are assumed to be complete before 
general shear begins. For the purpose of this analysis, general shear 
is defined as a sliding of a mass of soil along a shear surface, which 
occurs without a significant volume change. 

RESISTANCE DUE TO VOLUME CHANGE 

General 

For cases of track formation, where general shear takes place the 
maximum soil resistance available will be given by the resistance due to 
general shear plus resistance due to soil inertia; it is assumed that 
volume change and shearing distortion need not be considered. The greater 
part of the shearing distortion resulting from the rearrangement of 
particles and soil structure takes place prior to the onset of general 
shear. That is, when a state is reached at which the soil will fail by 
general shear rather than further shearing distortion, resistance to 
failure is no longer a function of shearing distortion. Except for very 
loose soils, a condition of general shear is forced in nearly all rolling 
sphere-soil slope interaction cases. On the basis of such reasoning, 

q v = 0- (5-8) 

Volume changes are, however, important in evaluating inertia effects 
of a shearing soil wedge during general shear; the looser the soil, the 
shorter will be the displacement along the shear plane, since more of 



the deformation associated with rolling has been taken up by shearing 
distortion and volume change. To be able to evaluate displacements along 
the shear plane at various soil densities , a relationship between 
sinkage and volume change is needed. 

Relationship Between Sinkage and Volume Change 

It is assumed that volume changes are complete before general shear 
begins. Figure 5-5 shows the geometry for the volume of soil displaced 
during densification and before general shear. 



Fig. 5-5. Geometry for analysis of soil displaced 
before general shear. 


The symbols in Figure 5-5 are defined as follows: z = track 

v 

depth or sinkage for maximum volume change before general shear 

z) / w y = track width for maximum volume change before general 

shear (w < w) . 
v 

The volume of soil displaced, Av , per unit length in the track 



direction is given by the segment of the circle below the original soil 
surface in Figure 5-5, and equals 


AV 


i 2 (^ - - 2e ) • 


(5-9) 


where 


0 = sin 

v 


“(t) • 


(5-10) 


Substituting Equation (5-10) into Equation (5-9) , Av can be expressed in 
terms of the track width to diameter ratio: 


* - 1 2 

7T . _ X / W . 

|^ Sln \1 

r) - sin 

2 “"“(t) 



l 


L J 

i 


(5-11) 


In terms of the sinkage, z^, it is also easy to show from Figure 5-5 that 


0 = cos 

v 


'(■ • *) ■ 


(5-12) 


Then the volume change as a function of sinkage becomes: 


A V = =r 


1 

r n 

1 

1 _ / z \ 

/ z \ 


— — — — roe’Hl — — 1 — c=in 



190° COS l 1 r j 

£ COS 1 1 

\ r / 




1 


(5-13) 


While it may not be possible to determine z^ exactly, it can be 
estimated from Equation (5-13) by assuming that 


A v = V - V , 
w r 


(5-14) 


where V = volume of track below original soil surface, and V = volume 
w r 

of rims above original soil surface. Equations (5-11) and (5-13) are 
plotted in Figure 5-6. The quantities and are measured from track 
cross sections. 

RESISTANCE DUE TO GENERAL SHEAR 
Selection of the Failure Surface 


In analyzing soil failure on the basis of a limit design method, 
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such as general shear, it is necessary to select a realistic failure 
surface. 


As shown in Chapter 4 , general shear appears to be related to the 
maximum pressure, which is about twice the average pressure along the 
sphere-soil contact. Therefore, the contact area controlling general 
shear is considered to be described by a smaller area than the total 
contact area. A contact area controlling general shear defined by r0/2 
in the longitudinal section as shown in Figure 5-7a was selected 
because it establishes a failure surface which agrees with experimental 
data (Figure 4-21a) . 

The sphere-soil contact or boulder-soil contact is usually rough. 
Meyerhof (1955) found that the base angle, ip, (see Figure 5-9a) for a 
rough strip foundation on cohesionless material was close to 1.2 <J>. A 

line drawn through distances r6/4 laid off along the sphere surface and 

/ 

vertically below, as shown in Figure 5-7b, defines an angle ABC close 
to 1.2 <f ) . This defines the beginning of the failure surface directly 
beneath the sphere. In the zone of radial shear, the failure surface 
in sand is a logarithmetic spiral (log-spiral) . This is a consequence 
of the requirement that slip lines intersect at an angle of 90° ± (p 
(Scott, 1963). In the passive Rankine zone, the failure surface is - 
a plane. 

In analyzing boulder track phenomena, the only measurable quantities 
usually are the track width, w, the boulder diameter, D, and the slope 
angle a. Having measured these quantities, the following additional 
steps establish the failure surface: 


1. Estimate an initial value for the angle of shearing 
resistance, <p. 


2. Calculate 0 and 0 (Sokolovski, 1960; Karafiath and Nowatzki, 1968) : 

2 3 


6 = 2 
2 2 


tt , . . -i sxn a 

— - ffl - a + sin — : t- 

2 Y sin 4> 


(5-15) 


e = 7 - <t> - e . 

32 2 


(5-16) 




It is to be noted that 0^ and 0^ are for the slope angle, a, 
without considering the bow wave. Experimental evidence on 
track formation (Chapter 4) shows that the bow wave is for many 
cases as high above the original surface as the bottom of the 
track is below. Figure 5-7 was drawn to reflect this, and it 
was assumed that 0^ and 0 are adequately representative also 
for the somewhat steeper slope angle caused by the bow wave . 

3. Extend the 0 line to intersect the bisector of the contact area 

2 

this establishes the center of the log-spiral. (The center of 
the log-spiral is assumed to be located as described because 
this location gives a realistic failure surface supported by 
experimental data discussed below.) 

4. Determine the log— spiral surface, and draw the failure surface 
to exit the slope at an angle 0^ as shown in Figure 5-7b. 

To see how well a longitudinal section of the failure surface, se- 
lected by the procedure outlined above, describes the true failure 
surface, predictions were compared with actual measurements. Figure 5-8 
shows the results of this comparison. Predicted distance of forward 
shear (the distance at which the failure surface exits the slope) is 
plotted in comparison with measured distance of forward shear (Figure 
4-21a) . The measured distance of forward shear is the distance from the 
crest of the track to the furthest shear surface; it is the width of the 
forward shear zone. The measured distance of side shear is also shown. 
It is the width of the shear zone on the side of the track measured 
from the crest of the track. 

The selected failure surface appears to be adequate. The failure 

surface as shown in Figure 5-7b also compares well in shape with those 

observed in model studies (Chapter 3) . The ratio of measured side shear 

to measured forward shear is close to 0.38 for the range of v /v ratios 

r s 

shown. This value is the experimental basis for the width of the 
failure zone. 


DENSITY RATIO 


O MEASURED SIDE SHEAR 
• MEASURED FORWARD SHEAR 
A PREDICTED FORWARD SHEAR 



RATIO OF MAXIMUM WIDTH OF SHEAR 
ZONE TO AVERAGE TRACK WIDTH 


Fig. 5-8. Density ratio vs predicted and measured width of 
shear zone. 




95 


Resistance Due to General Shear 

Previously (Hovland and Mitchell, 1970) a method was developed for 
evaluating boulder track phenomena using modified bearing capacity 
theory (see Chapter 2). The general bearing capacity equation, (2-4), 
was first derived for an infinitely long footing on a horizontal surface 
by Terzaghi (1943) . The derivation was accomplished by solving separately 
for the contributions of density, cohesion, and surcharge. These con- 
tributions were added to get the total soil resistance by assuming 
superposition to be valid. It has since been shown that superposition 
is not exactly valid for this case (Scott, 1963) , but that it is a very 
good assumption and that Equation (2-4) is in close agreement with 
experimental and field data. 

Subsequently the bearing capacity factors (N , N , N ) in Equation 

I '■i 

(2-4) have also been determined for an infinitely long footing on a 
slope (Meyerhof, 1951). 

The general bearing capacity equation has been modified to make 

it applicable to footings of various shapes other than an infinitely 

long footing. The approach has been to adjust the basic equation, (2-4) , 

by shape factors (s , s , s ) as described in Chapter 2. It has been 
J y c q 

found from experimentation that a change in the shape of the footing is 
much less significant than a change in <J> on the resulting bearing 
capacity. Adjusting Equation (2-4) by shape factors is believed to be 
logical and generally valid. 

In solving the rolling sphere-soil slope interaction problem, the 
actual semicircular contact area was first approximated with a rectan- 
gular area of equal size (Chapter 2). By using shape factors. 

Equation (2-4) was modified for an approximation of the rolling sphere- 
soil slope interaction problem. The resulting expression was 

^ . 0.188 N y +1.1 + 0.55 (5-17) 

In Equation (5-17) q g is the unit bearing capacity in earth gravity and 
y g is the soil unit weight in earth gravity. N^, N^, and are the 
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bearing capacity factors. 

Equation (5-17) predicts soil resistance due to general shear for 
the sphere-soil contact area as if the failure surface was as predicted 
from bearing capacity theory for a footing. Based on experimental 
evidence (Chapters 3 and 4) , a method of predicting a more realistic 
failure surface for the rolling sphere-soil slope interaction problem 
was described in the previous section. It is now possible to modify 
Equation (5-17) to reflect the more realistic failure surface. 

Cohesion acts only along the failure surface. Therefore, for the 
same conditions and same bearing area, the ratio of the failure surface 
area for the rolling sphere, S , to that of a footing, S, , is an index 
of the extent by which the cohesion contribution to resistance to 
rolling differs from the contribution of cohesion to bearing capacity 
under a footing. A footing resting on a horizontal surface with a 
usual failure surface is illustrated in Figure 5-9a. The failure sur- 
faces for a rolling sphere and a footing for analogous conditions are 
illustrated in Figure 5-9b. (S s and were measured as the traces of 

the failure surfaces in the vertical plane.) The ratio S /S was deter- 

s f 

mined graphically for several conditions of w/D, <(>, and slope angle a. 

The results are plotted in Figure 5-10. In the comparison, a ^ = 45 + <f>/2 
was used for the footing; this value of \p was used by Meyerhof (1951) in 
determining the bearing capacity factors. 

From Figure 5-10, 

S 

= 0.37 + 0.25 . (5-18) 

► 

Therefore, from the proportionality. 


where N 
sphere. 


cs 


N S 


cs s 



N = f 0.37 + 0.25 , 

cs V D/c 

is the bearing capacity factor for cohesion for 


a rolling 


(5-19) 

(5-20) 


w 



a) Failure surface for a footing as indicated by bearing 
capacity theory. 



b) Failure surfaces for a sphere and a footing for analogous 
conditions as indicated by bearing capacity theory. (For 
failure condition on a plane normal to the direction of 
travel, see Figure 5-13). 


Fig. 5-9. Illustration of the failure surfaces for a rolling 
sphere and a footing. 
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Surcharge acts on an area confined by the failure surface and is , 

therefore, a function of the extent of the failure surface similar to 

cohesion. The relationship between N and N has been established and 

c q 

can be found in many books and papers on the subject (Terzaghi, 1943): 


N c = cot 4> - lj . (5-21) 

Therefore , 

N qs * (o.37 + 0.25 jA tan <J> N, + 1, (5-22) 


where N^ g is the bearing capacity factor for surcharge for a rolling 
sphere . 

Resistance caused by the weight of soil, the term, is a 

function of the longitudinal section area of the failure wedge. The 

ratio of the longitudinal section area of the failure wedge for a 

rolling sphere, A , to that of a footing, A , may therefore be taken as 

an index of the extent by which N differs from N . (N is the 

ys y ys 

bearing capacity factor for soil weight for a rolling sphere.) These 

areas are illustrated in Figure 5-9b. The ratio of A /A, was deter- 

s f 

mined graphically for the same conditions of w/D, <J>, and a used to 
determine S s / S £* The results are plotted in Figure 5-11. A comparison 
of values of S g /S f and A g /A f in Figures 5-10 and 5-11 show that 



Therefore, from the proportionality 


N A 
Ys _ _s_ 

\ "*£ 


N ys - (°' 37 + °‘ 25 d) 2 \ ' 


(5-24) 


(5-25) 
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Replacing the bearing capacity factors in Equation (5-17) with 

N , N , and N , we have 
ys cs qs 

- 0.188 N I + 1.1 N +0.55f-'\w . (5-2 

W Y S ys \wy s f CS J qs 

Equation (5-26) , therefore, provides an expression for rolling sphere- 
soil slope interaction based on the selected failure surface. 


Summary 

In summary , the steps and assumptions leading to Equation (5-26) 
are reviewed. 

1. The general bearing capacity equation (2-4) which was derived 
by Terzaghi (1943) , is the basis for Equation (5-26) . Equa- 
tion (2-4) was developed for an infinitely long footing, which 
implies a two-dimensional analysis. Superposition of the 

effects of N , N , and N was also assumed, 
y c q 

2. Equation (5-17) was developed from Equation (2-4) using shape 
factors . 

3. Studies of the failure mechanism (Chapter 3) and the pressure 
distribution (Chapter 4) showed that the actual failure wedge 
(zone of shearing soil) in the rolling sphere problem is much 
smaller than that predicted for a footing of equal area from 
bearing capacity theory. It was further shown that shear 
failure is more directly related to the maximum pressure than 
to the average pressure underneath the sphere. These two 
characteristics were considered in this section, and on that 
basis. Equation (5-26) was developed from Equation (5-17). 

Therefore, Equation (5-26) was developed from a two-dimensional 
analysis and modified to be applicable to the rolling sphere-soil slope 
interaction problem. It is believed that it closely approximates the 
actual three-dimensional case. Theoretical predictions are compared 
with experimental data in the next chapter. 


102 


RESISTANCE DOE TO SOIL INERTIA 

It is possible to speak of two types of motion for a rolling 
sphere: uniform rolling and bouncing. Apparently the effects of even 

minor irregularities in the soil slope and sphere geometry become mag- 
nified as the velocity increases. At a certain velocity, depending on 
the situation, conditions of uniform rolling are influenced to the 
point that bouncing begins. The theoretical developments in this chap- 
ter are only for uniform rolling. 

Experimental investigations on accelerating spheres have shown 
that the track width is essentially independent of sphere velocity until 
the sphere starts bouncing. Therefore, one might suspect that inertia 
effects of the moving soil would be insignificant. On the other hand, 
because of the reduced lunar gravity, inertia effects on the moon will 
be six times more important, relative to soil resistance due to fric- 
tion and surcharge, than on the earth. Thus, inertia forces must be 
evaluated in order to determine their importance. 

The resistance due to the inertia force of the moving soil, F , is 

the product of the mass of soil involved, m , and its acceleration, a ; 

s s 

F = m a . (5-27) 

I s s 

These items will be analyzed below. 

Mass of the Moving Soil 

Due to mathematical complexities in expressing exactly the volume 
of soil bounded by the log-spiral failure surface, an approximation of 
the volume of soil is developed. Should the volume of soil turn out 
to be significant, a more precise determination may be necessary. 

A longitudinal section of the failure surface confining the moving 
soil, and an approximation to that surface confining approximately the 
same area of soil as shown in Figure 5-12a. The ratio of d/d is 
plotted in Figure 5-12b; this ratio was determined graphically. 
Mathematically, the curve in Figure 5-12b can be expressed by 
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a) Boundary confining approximately the same area of 
soil as the actual failure surface. 



b) Ratio of d/d^ as a function of d>. 

o T 


Fig. 5-12. Geometry for volume of moving soil. 
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f” 3 1 + 0.002<t> 1 - 77 , (5-28) 

o 


where <{> is in degrees. 

From Figure 5-12a, the volume of soil is approximately 


V 

s 


2 tan 0 


* 


(5-29) 


where b is the width of the moving soil wedge, 
s 

From Figure 5-8, the maximum width of the side shear zone is 

shown to be 0.45 w for y /y =1.7. The depth of the failure surface 

below the bottom of the track can be estimated graphically to be about 

0.6 d Q from drawings similar to that shown in Figure 5-12a for various 

w/D ratios. With these dimensions and a w/D ratio of 0.6, a cross 

section for the maximum depth and width of the moving soil wedge can be 

drawn as shown in Figure 5-13. As shown in Figure 5-13, the width of 

the failure wedge for purposes of computing can be approximated with 

b = 1.4 w. 
s 

From Figure 5-12a 


d 

o 




r 

k sin 

2(2 V 


w 

D * 


(5-30) 


Substituting into Equation (5-29) for b , d, and d , and dividing 

so 

through by r w results in a dimensionless expression for the volume of 
moving soil: 

V s 0 ^7 (l + 0.0Q2(t> 1 - 77 ) / . -i w \ 2 

r 2 w " 8 tan 0 3 \ Sin D / ’ (5-31) 


Equation (5-31) is plotted in Figure 5-14. The approximate mass 
of moving soil is then found simply by multiplying the volume by the 
soil unit weight and dividing by the acceleration of gravity: 

s_ 

g 


m = 
s 


(5-32) 


5 
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Fig. 5-14. Volume of sheared soil wedge as a function of 
w/D and 4> . 
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Acceleration of the Moving Soil 

Relationships between the velocity of the sphere and the velocity of 
moving soil; 

The acceleration and the velocity of the moving soil must be 
related to the velocity of the sphere. Since the velocity of the sphere 
can be readily determined experimentally, it is desirable to express 
the velocity of the moving soil in terms of the velocity of the sphere. 

Referring to Figure 5-13, continuity requires that 

Vt = vA w ' (5 “ 33) 

where v g is the velocity of moving soil, v is the velocity of the 
sphere, is the total cross-sectional area of moving soil, and A^ is 
the relevant cross-sectional area of the track. 

For loose sand, compression of the soil beneath the sphere reduces 
the magnitude of displacement along a failure surface. The relevant 
track area is, therefore, somewhat smaller than the track area associ- 
ated with the crest- to-crest track width. Both the track area 
associated with the crest-to-crest track width and the reduction due 
to volume change (compression) can be expressed in the form of Equation 

(5-9) to give for A , 

w 

A = — 0 - sin 2©\ - (— 0 - sin 20 ^ * (5-34) 

w 8 \90 J 8 \90 v w) 

or 

(5-35) 

The total cross-sectional area of the moving soil, A fc , can be 
expressed from Figure 5-13 as 

A. = 1.4 w (0.6d + z) . (5-36) 

t o 

Combining Equations (5-33) , (5-35) , and (5-36) , the velocity of 


A = (0 - 0 \ - sin 20 + sin 20 1 

w 8 L90 \ v/ vj 
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of the moving soil is 

= \r v = I ( 9 ' e v) ' Sin 26 + Sin 26 y] 

1.4w (o.6d Q + z^ 


v = — v 
S A t 


v , 


(5-37) 


since 


j ~ / o\ ^ 

d O - — (2) ' 


ft . -1 w 

0 = sin - , 


w 

6 V = sin -1 , and 


z = r^l - cos 0^ = r. 1 - 1 - 


cos QJ = 

Equation (5-37) can also be expressed as 


„ / «\ / \ / » \ 

IT I . _i W . -1 V | I ft • -1 W I . I „ . -1 V I 

— { sin — - sin —I - sin \ 2 sm — I + sin 1 2 sin — / 

? 90 \ D D / \ D/ V, D / 

In order to appreciate the relative magnitudes of and v. 


v . 


(5-38) 


Equation (5-38) was solved for w/D = 0.69 (average value for the lunar 
boulder tracks analyzed) , and assuming that Av = 0 (meaning that 
w^ = 0) . This gave 


v = 0.32 v . (5-39) 

s 

( An alternative approach for estimating the velocity of the moving 

soil results from consideration of a longitudinal section of the 
rolling sphere-soil slope problem. Figure 5-15 shows the conditions 
to be analyzed. 
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If 



Fig. 5-15. Sphere and soil movements. 




In the time required for the sphere to move from A to B, t S p here 
the soil at C must have moved from C to D. (That the soil moves 
approximately from C to D was demonstrated experimentally from model 
studies , see Figure 3-7.) Due to compression of the soil (volume 
change) , the displacement generated along the failure surface will be 
somewhat smaller than the distance CD. This displacement can be ex- 
pressed by (z — z ). Equating the time intervals, 

( z - 2 v) 


t sphere ^oil 


9r 

v 


(5-40) 


from which 


■ (‘ - z v) 57 


(5-41) 


Equation (5-41) expresses a relationship between the magnitudes of v g 
and v; it says nothing about their directions, which are irrelevant in 
this consideration. 
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Since 0 = sin -1 jj- and z = r * 




we also have 



(5-42) 


In order to compare Equation (5-42) with Equation (5-38 ) , 

Equation (5-42) was solved for w/D = 0.69 and z^ = 0 (meaning w^ = 0) . 
This gave 

v. = 0.36 v. (5-42b) 

• 

Since Equations (5-38) and (5-42) give .nearly the same results as 
shown by Equations (5-39) and (5-42b) , and since both are approximate, 
Equation (5-42) will be used because of the simpler calculations 
involved . 

Equations (5-41) and (5-42) give the velocity of the shearing soil 

wedge as a function of the velocity of the sphere. The and z terms 

can be determined from Figure 5-6 which applies to both, although 

marked out for z . 

v 

Expressing the distance through which the soil moves by (z - z^) 
as done above, is an approximation and implies three situations: 

1. When there are no raised rims; Av > 0, = z, and there is 

theoretically no movement of a mass of soil since there is no 
general shear. This case is encountered only in very loose 
soils. 

2. When the volume of the rims equals the volume of the track; 

AV * 0, z = 0, and the soil moves approximately through the 
distance (z - z^) . 

3. When the volume of the rims is greater than the volume of the 
track (a dilating soil) ; assume that Av - 0, and z^ = 0 as 
above for case 2. 
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Soil acceleration: 

Since the soil wedge moving during general shear starts from rest 
and comes to rest after the rolling sphere has passed , we know that the 
initial and final soil velocities are zero. Somewhere in between, the 
velocity of the shearing soil reaches a maximum. As an illustration, 
for a parabolic velocity function, the displacement, velocity, and 
acceleration as a function of time would plot as shown in Figure 5-16. 



Fig. 5-16. Time functions of x„. v„, and a.. 
* s s s 


If the velocity of the soil reached a maximum before t /2, the 

SO 

acceleration function would be curved as shown by the dashed line in 
Figure 5-16. From a mathematical point of view, it has been found most 
convenient to represent the acceleration curve by 


a g = -Bjt 1 ^ 11 + B z , where Bj and B 2 are 


(5-43) 


constants to be determined from the boundary conditions, and n is a 
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positive number. Therefore, at t « 0, a = a. , initial acceleration, 

S X 

giving B = a.. At t = t , a = a _ , final acceleration (actually 
21 so s f J 

deceleration) giving 


a = - B t ly/n + a. , 
f l so x 


B . t l/n 

SO 


(5-44) 


Then the acceleration becomes 


■ ai ' ( a ‘ ' ^ (t)‘ 


(5-45) 


The velocity is found by integrating Equation (5-45) 


v ■ f a dt = a . t - / a . - -—2 — - — 

s Is 1 \ 1 f/ 1 + n l 


(i+n)/n 

~J7n 


+ B , (5-46) 

3 


where B is an integration constant. At t = 0, v =0 giving B =0. 

3 S3 

Also, att=t ,v =0 giving 
so s 3 ” 


0 = a.t - (a. - a,') t — — t , and 

iso \i f/l + n so 


(5-47) 


since t / 0 , 
so 


a . a . 

a ^ or n = — 

f n a f 


(5-48) 


After some manipulation of Equation (5-46) , 


r it v /n i 

•v i -(H 

L \ so/ 


\ (5-49) 


The displacement can be found by integrating Equation (5-49) 




+ B , where (5-50) 

4 
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is an integration constant. At t = 0, x g = 0 giving = 0. Then 


■ ai * 2 D ‘ 1 ” 2» (tj ]• 


We have yet one more conditions at t = B t , v = v , and 

5 so s max 

a =0. Hence , 
s 


[■ - Hr^ r;) ]• 


since a. ^ 0 
x 


Then: 


, a nd 


■.-w- 

v ■ n (ittt) So 1 - [(rrr) 

, ( . ,\ /n + l\ n V max 

+ ) — • 


Substituting a^ from Equation (5-55) into Equations (5-45) , 

(5-49) , and (5-51) results in the final general expressions for accel- 
eration, velocity, and displacement respectively: 

\<« - ” {Hrf" - Hr (F)' /n ] (5 ' 56 ’ 

- (" ♦ *) M t [ - feT] <5 - 571 


* s (t) 


■ (n + x) (^F ' ^(t:) ,/n ] v ’ 


We are, however, interested in the total or maximum displacement. 
This can be obtained by integrating Equation (5-57) between known 



114 


limits (0 and t ) . Then the total displacement equals 

SO 


total 


- n + 1 / n + l\ n 
2(1 + 2n) \ n / V ] 


t 

max so 


(5-59) 


The max i mum inertia force develops for maximum acceleration. It 
was shown previously that for the function assumed, a g = a max when 
time equals zero. Therefore, from Equation (5-56) 


max 


=(— )(^n 


max 

so 


(5-60) 


Since we are only able to infer the total displacement from mea- 
surements of volume change and track dimensions , which implies a 
relationship with the average velocity, a theoretical expression for 

a in terms of average velocity is needed. The theoretical average 
max 

velocity is 


V at 


L 


so 


v dt 
s 


so 


X max _ n + 1 / n + l\ n v 

t 2(1 + 2n) \ n / max ’ 

so 


(5-61) 


This velocity, v , can be related to the velocity expressed by 
at 

Equation (5-41) , v g : 


v = v . 
s at 


/ \ v n + 1 / n+l\ n 

V Z v) 9r 2(1 + 2n) \ n J V max* 


(5-62) 


Solving for v , 
’ max 


( z ~ z v) 


9r 


max 


n + 1 /n + l\ n 
1 + 2n) \ n / 


2(1 


(5-63) 


Dividincr v by t = 9r/v, as expressed by Equation (5-40) , 
3 max ■'so V 


v max Z 2 v / v \ 

So n>l 


(5-64) 


2(1 + 2n) 
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Substituting Equation (5-64) into Equation (5-60) results in the 
desired expression for maximum acceleration: 

%ax - 2(1 * 2n> lz - v(e7) 2 ' (5 - 65) 

where again 0 = sin” 1 . 

Equation (5-65) is a simple expression where the terms in the last 
two parentheses can be experimentally determined. A suitable means for 
determining n has not yet been developed? however , some experimental 
data which are useful for its estimation are presented and discussed 
in the next chapter. Since n = ~ a ^/ a f / Equation (5-65) also states 
that the closer the maximum velocity is to the start of motion timewise, 
the higher is a * 


Inertia Force Associated with the Moving Soil 

Combining Equations (5-32) and (5-65) , we have for the inertia 
force: 


F 


I 


V y 

2(1 + 2n) 
9 



(5-66) 


This force can be expressed as force per unit area analogous to unit 
bearing capacity by dividing by the contact area A = (tt/2) (w/2) 2 : 


q_ = ~ (1 + 2n) 

I 7T 



(5-67) 


A convenient dimensionless expression results from rearranging the 
terms in Equation (5-67) : 


V z ' V 


V 

s 




(5-68) 


To show that the inertia resistance is independent of gravity, 
Y can be expressed as pg, where p is the mass density of the soil. 

S 


Then Equation (5-68) becomes 


p(z - z v ) 

I 

In solving Equation (5-68) or (5-69) , V is determined from 

s 

Figure 5-14, and z and z^ can be determined from Figure 5-6. As an aid 
in solving Equation (5-68) it is plotted in Figure 5-17. Comparison 
of relative magnitudes of qj and q g are presented in Chapter 7. For 
the solid 12-cm spheres rolled at WES, q as a percentage of (q + q ) 
varied from about 1% to 16% for moist Yuma sand and dry Yuma sand, 
respectively. 


w 2 D 


= 20.3 


(1 + 2n) 
!) 2 


D 


(5-69) 


SUMMARY 


It is now possible to combine the effects of general shear and 
soil inertia into one equation; this presumes that superposition is 
valid. First, it is necessary to rearrange Equation (5-68) to obtain 
a term comparable to the terms in the bearing equation. Such 
rearrangement gives 




B (z - z ) 

v 


V 


w 3 D 


t 


(5-70) 


where 


B = 20.3 


(1 + 2n) 




(5-71) 


as obtained from Figure 5-17. 

Then, the total unit pressure resistance to rolling is given by 



0.188N + 

ys 




(5-72) 


In Equation (5-72) q g is the unit resistance in earth gravity. 
In lunar gravity, replacing the y g term by Y g /6, the unit resistance 
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Fig. 5-17. Plot of Equation (5-69). 



is given by 


wy 


■« 0.0314N + 1.1 

ys 




N + 
cs 


0.0916 


( 5 )'. 


qs + wy 


(5-73) 


From the relationship of sphere weight to bearing area (Hovland 

and Mitchell, 1970), a second expression is obtained. This expresses 

/ 

the average applied pressure, and in earth gravity is 


%_ 4 (v0 
"*8 ‘ 3 (»/.>)’ ' 

In lunar gravity 

%5_ 2 (vO 

“*8 ' 9 (w/d) j ’ 


(5-74) 


(5-75) 


Further, as shown in Figure 5-9b, a unique relationship must exist 
between w/D and the slope angle, a, for constant velocity rolling. 

The relationship is 

jj = sin 2a . (5-76) 

We, therefore, have three equations which can be combined to solve 
for parameters of interest. These equations are: 

1. Equations (5-72) or (5-73) expressing total soil resistance, 

2. Equations (5-74) or (5-75) expressing the applied pressure, 
and 

3. Equation (5-76) ejqpressing the relationship between w/D 
and a. 

The parameters associated with these equations , known or unknown 
depending on the situation, are: 


1. y , y , <J>, c, and v (usually variables) 
s r 



2. w, D, Av, a, n, and 9 (usually measured, known, or assumed) 

3. z, z , V ,q , and q (usually calculated from developed 

V S G X 

relationships) . 

Graphs for easy solution of the above equations can be prepared. 
It is necessary first, however, to examine experimental data in terms 
of the theory. 



Chapter 6. ROLLING SPHERE EXPERIMENTS AND COMPARISON . 
WITH THEORETICAL PREDICTIONS 


In order to study simultaneously the influence and interrelation- 
ships of a number of pertinent parameters, solid spheres were rolled on 
Yuma sand at the U. S. Army Engineer Waterways Experiment Station (WES) 
in Vicksburg, Mississippi. Because of the extensive facilities and 
excellent cooperation provided by the WES personnel, it was possible for 
us to extend these experimental studies beyond what was originally 
planned. The assistance and facilities provided by WES to this effort 
were funded through NASA Defense Purchase Request H-58504A from MSFC. 

This chapter presents a description of the equipment used , the 
testing technique, and the results of these experimental studies. The 
results are presented in summary form; additional details and data are 
to be found in the appendices. 

The analysis of the experimental data can be divided into 1) measure- 
ments of the tracks left by the rolling spheres, and 2) analysis of films 
taken of the rolling spheres. 

Wherever possible, curves developed according to the theory in the 
previous chapter are superimposed on the experimental data, with a brief 
discussion on the significance of the comparison. 

DESCRIPTION OF EQUIPMENT AND TESTING TECHNIQUE 

The testing was done on Yuma sand prepared to uniform density and 
moisture content in a large soil car. Methods for preparing the soil 
car, controlling moisture content, and determining soil properties are 
described in Appendices A and B. The cars were 27 feet long, 64 inches 
wide, and contained soil to depths of 18 or 32 inches. Figure 6-1 is a 
photograph of the car used. 
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Fig. 6-1. Soil car used for testing at Vicksburg. 





Nine soil cars were tested, with soil conditions as indicated in 
Table 6-1. 


Table 6-1. Sand properties in different soil cars. 


Car 

# 



0) 

% 



G 

(MN/m 3 ) 

G 

(lb/in 3 ) 

^d 

(gm/cc) 

<P 

(°) 

c 

KN/m 2 

c 

lb/in 2 

D 

(%) 

1 

0.20 

0.57 

2.11 

1.480 

37 

0 

0 

34 

2 

0.20 

0.45 

1.67 

1.453 

37 

0 

0 

26 

3 

0.20 

3.02 

11.20 

1.618 

42 

0.20 

0.029 

86 

4 

1.50 

2.65 

9.82 

1.473 

37.5 

1.07 

0.156 

35 

5 

1.55 

0.20 

0.74 

1.245 

32 

0.46 

0.067 

< 0* 

6 

1.55 

0.24 

0.89 

1.245 

32 

0.46 

0.067 

< 0* 

7 

1.50 

0.35 

1.30 

1.280 

33 

0.47 

0.069 

< 0* 

8 

1.50 

1.90 

7.04 

1.491 

37.5 

0.84 

0.122 

41 

9 

1.55 

1.76 

6.52 

1.436 

36.5 

0.88 

0.128 

20 


In Table 6-1, the symbols at the columns are defined as follows: 
w% = moisture content in percent 

G — cone index gradient expressed in mega newtons per meter 
cubed and in pounds per cubic inch 
Y d = dry density of sand in grams per cubic centimeter 
<t> = angle of shearing resistance (friction angle) , determined 

from Figure B-10 and based on triaxial results 
c = apparent cohesion expressed in kilo newtons per meter 
squared and in pounds per square inch 
D = relative density of sand in percent 


* D r = 0 refers to the loosest state of air-dry Yuma sand; moist Yuma 

sand can be prepared to an even looser state with respect to Y , 

d 

therefore giving D r < 0. 



Eight different spheres were rolled over the soil in cars 1 through 
5, and selected spheres were tested with cars 6, 7, 8, and 9. Two series 
of data were collected: 1) track formation from spheres of varying 

density but constant size, and 2) track formation from spheres of 
varying size but constant density. Data from track measurements are 
presented in Appendix D. 

The spheres were made of concrete, giving a rough contact between 
soil and sphere. They were cast in hemispherical molds and subsequently 
the two halves were cemented together using an epoxy cement. The 
desired variations in density were obtained by using lead shots as a 
heavy aggregate and styrofoam plastic as a light aggregate. Sphere E, 
the lightest sphere, was made from a piece of extremely light weight 
concrete. Characteristics of the spheres are listed in Table 6-2. 


Table 6-2. Sphere characteristics. 


Series 

Sphere 

Diameter 

(cm) 

Density 

(gm/cc) 

Density, 

A 

12.10 

3.69 

(constant size) 

B 

12.15 

2.14 


C 

12.00 

1.73 


D 

12.10 

1.05 


E 

12.02 

0.52 

Size 

la 

7.16 

1.86 

(constant density) 

2(c) 

12.00 

1.73 


3a 

17.40 

1.76 


4 

25.30 

1.79 


Movie films were taken of all rolling spheres using a Bell & 

Howell camera with 1", 2", and 3” lenses. The camera also had variable 
speed adjustment, and we used speeds of 25 and 50 frames per second 
depending on the focus and lighting conditions. The movies provided 
information on the distance, velocity, and acceleration of the spheres 
as a function of time. Equally significant was the information the 
movie films provided on the rolling sphere — soil slope interaction 
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phenomenon and on actual movement of the shearing soil wedge. The film 
is on file in the University of California, Geotechnical Engineering 
Laboratories. The films were studied using a film analyser which allowed 
observation frame by frame. The positions of the camera and soil car, 
for both cars tilted sideways and cars tilted on edge, are shown in 
Figures 6-2 and 6-3. With this geometry , distance to the sphere and the 
scale of any frame could be determined if necessary. 

Figure 6-1 shows the soil car tilted to the desired slope, the 
camera, and the sphere in position retained by a "flip gate". The flip 
gate was specially designed to open the instant the camera was started. 
This was achieved by a solenoid on top of the flip gate. The solenoid, 
the moment triggered, would release the gate allowing the sphere to roll. 
A detailed drawing of the flip gate is shown in Figure 6-4. 


TRACK MEASUREMENTS 

Slope Angle Required for Rolling 

Records were kept during the testing of the slope angle at which a 
sphere of a particular density started rolling. This relationship was, 
of course, a function of the soil conditions. The results are presented 
in Figure 6-5. 

In the previous chapter. Equation (5-76) was developed for track 
size as a function of slope angle for constant velocity rolling: 

jj = sin 2a, (5-76) 

The theoretical curve for this Equation is shown in Figure 6-6. In order 
to plot the experimental data on the same figure (Figure 6-6) , the 
straight lines on Figure 6-5, separating rolling and no rolling situations, 
were combined with experimental information on the relationship between 
the density ratio, y^/ y^, and the w/D ratio. Experimental data for the 
relationship between Y r /Y s and w/D are presented in Figure 6-7. 

Therefore, the experimental verification of Equation (5-76) was 
achieved for certain soil conditions, by first finding the y /y ratio 

IT S 

for a selected slope angle from Figure 6-5 and then finding the w/D ratio 


CENTER OF CAMERA LENS CRANE 
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Top view 
Scale: 1" = 5" 


Fig. 6-4. Detail of flip gate 






DENSITY RATIO r / v DENSITY RATIO 


28 



0 10 ° 20 ° 30 ° 40 ° 

a) Slope angle, a, required for rolling in dry sand 



0 10 ° 20 ° 30 ° 40 ° 

b) Slope angle, a, required for rolling in moist sand 


Fig. 6-5. Density ratio vs slope angle required for rolling on Yuma sand 



L29 


Slope angle limitations: 

1 . a < <}> 

2. For constant velocity rolling a and w/D are 
related by: 

w/D = sin 2a, or a = 0/2 


THEORETICAL CURVE 
-jj- * SIN 2 a 


LEGEND 

• CAR 4, MOIST DENSE SAND *<*J = 1.50 % 
O CAR 3 , DRY DENSE SAND 
+ CAR 5 , MOIST LOOSE SAND ,0) = 1.55% 
A CARS 18 2, DRY LOOSE SAND 


Fig. 6-6. Slope angle, a, required for rolling as a function of w/D. 
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for that Y r /Y g ratio for comparable soil conditions from Figure 6-7. It 
is to be noted that Figure 6-7 contains experimental data for many slope 
angles; therefore, the comparable condition is found by matching not only 
the y/y ratio but also the slope angle, a, on Figure 6-7 as closely as 
possible with the slope angle selected at the start on the basis of the 
trend of data from Figure 6-5. Experimental curves were replotted on 
Figure 6-6 by this procedure. The agreement confirms the proposed rela- 
tionship, Equation (5-68). The agreement is best for loose moist soil. 

In this material the tracks left by the rolling sphere had sharp, well 
defined edges which could be accurately measured. The deviation is 
greatest for loose dry sand where the tracks tended to cave or slump and 
the crest-to-crest track width, w, measured was undoubtedly somewhat 
greater than the track width at the time the sphere was in contact with 
the soil. 

The comparison in Figure 6-6 supports the concept that the contact 

IT /w \ 2 

area between sphere and soil is equal to — y —I . 

Track Width as a Function of the Density Ratio 

The experimental data on the w/D ratio as a function of the Y r /Y g 
ratio and other soil conditions are presented in Figure 6-7. The crucial 
question was whether this experimental evidence, involving various soil 
(sand) conditions, could have been predicted from the theory. The 
appropriate equations, e.g. 



— = sin 2a 


(5-76) 



DENSITY RATIO , y '/ y . DENSITY RATIO 


3 


LECEND (ALL CARS) 

• O s | 0 « 

O « * 15° 

+ « * 20° 

AO* 25° 

X O * 3 5 * 


j > s 370 

Y» * 1-465 


+ O /t 


c « 0.2 KN/m 


a) CARS I AND 2 

( LOOSE AIR-DRY SAND) 


b) CAR 3 

( DENSE AIR- DRY SAND) 


9 • 37.5* 

Y% * 1.47 
C s 1.07 KN/m 1 

a 1 = 1.50% 


/ 


Y, » 1.23 
C * 0.46 KN/m 1 
to = 1.55% 


C) CAR 4 

(DENSE MOIST SAND) 


d) CAR 5 

(LOOSE MOIST SAND) 


Fig. 6-7. Density ratio vs track width over diameter ratio for spheres 
rolling on Yuma sand (curves are theoretical for a. required 
for rolling, points are experimental). min 
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were solved for the same soil conditions. The velocity term, v 2 /gD, was 
determined from data presented subsequently in Figures 6-12 through 6-15. 
The contribution due to soil inertia, was determined as outlined 

in Chapter 7. The results for smallest and largest values of w/D con- 
sidered, are presented in Table 6-3. 


Table 6-3. Relative contribution of soil inertia. 


Soil cars 

4> 

w 

*1 

Suggested reason for 

D 

x 1UU 

% + % 

contribution of 

If 2 

37° 

0.5 

16% 

Rel. high because c = 0 



0.9 

16% 

w 

and — was large 

3 

42° 

0.3 

12% 

Rel. high because cN was 

cs 



0.7 

20% 

small and velocity was fairly 
high 

4 

37.5° 

0.1 

1% 

Low because cN was high 

w ^ ^ 



0.5 

6% 

and ~ was small 

5 

32° 

0.4 

0.6% 

Low because soil was loose 



0.9 

1.3% 

and there was little or no 





general shear 


The theoretical curves determined as explained above are plotted in 
Figure 6-7. 

The agreement supports Equations (5-72) , (5-74) and (5-76) . It is 
to be noted, however, that the agreement may be coincidental to some 
extent since Equation (5-72) was not rigorously developed but contains a 
number of assumptions and approximations. Therefore, while the comparison 
in Figure 6-7 should be considered evidence for the validity of the 
proposed theory for the range of soil conditions and sphere sizes compared, 
significantly different conditions would at least require some empirical 
verification. The bearing capacity equation for footings, from which 
Equation (5-72) came, has extensive credibility in practice for its 
general validity. Therefore, having matched Equation (5-72) for a limited 
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number of sphere and soil conditions, it is believed that Equation (5-72) 
will properly predict many other conditions as well. 

Track Width as a Function of Sphere Diameter 

During experimentation it was noted that the w/D ratio was dependent 
on the size of the sphere. Theoretical predictions were made for this 
relationship using Equations (5-72), (5-74), and (5-76). The comparison 
between theoretical predictions and experiments is presented in Figure 6-8. 
The agreement is believed to be good. 

Volume Change as a Function of Sphere and Soil Densities 

Cross sections were measured for all tracks to investigate soil 
volume changes as a function of the density of the sphere and initial 
soil density. The cross sections were usually determined at a distance 
of 1.5 feet from start of rolling. Examples of track cross sections are 
to be found in Chapter 4, Figures 4—6 through 4—14. The average volume 
change was determined using Equation (5-14):* 

AV " V w " V r ' (5-14) 

where = volume of track below original soil surface, 

and V r = volume of rims above original soil surface. 

The sum of V and V was defined as V (V = V + V ) . 

w r t t w r 

The results of the track cross section measurements are presented in 

Figure 6-9. In this figure, volume increase (dilation) is positive, and 

volume decrease (compression) is negative. More data were available for 

a slope angle, a, of 20° or 25° than for other slope angles; therefore, 

lines were drawn through points for these slope angles as shown in 

Figure 6-9. 

The information in Figure 6-9 can be replotted in more familiar form. 
Since the lines in Figure 6-9 represent spheres of known densities, the 
AV/V ratio can be determined for assumed soil densities. This is done by 

1. assuming a soil density, y , 

2. computing the y /y ratio knowing y from Table 6-2, 

IT S IT 

* This method does not account for local volume changes within the regions 
considered. 


SPHERE DIAMETER , D , cm 


134 



Fig. 6-8. 


Track width over diameter ratio as a function of sphere diameter. 
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3. determining the resulting AVA fc from Figure 6-9 for the particular 
Y r (sphere) line, and 

4. plotting AV/V^_ vs. for the assumed soil density, as shown in 
Figure 6-10. 


FILM ANALYSIS 

Distance, Velocity , and Acceleration vs Time 

A measuring tape was placed alongside the anticipated track before 
the sphere was rolled. Knowing the speed at which the films, were taken, 
it was possible, with the aid of a film analyzer, to determine the distance 
rolled at any particular time. 

From markings on the sphere surface it was also possible to determine 
the number of revolutions rolled as a function of time. Slip could be 
determined as the difference between distance revolved and distance 
travelled (distance revolved as if the sphere had been rolling on a hard 
surface) divided by the distance travelled. 

The data were plotted as functions of time as shown in Figure 6-11. 

The velocity of the rolling sphere was determined as the slope of the 
distance vs. time curve. The acceleration of the rolling sphere was 
determined as the slope of the velocity vs. time curve. The data were 
adequate to reliably define the velocity of most of the spheres. Based 
on intermittant checks of the camera speed during testing, the velocities 
should be accurate to within 8% of the actual value. 

The remainder of the data giving distance, velocity and acceleration 
is on file at the University of California, Geotechnical Engineering 
Laboratory . 

Velocity and Acceleration vs Density Ratio 

Velocity and acceleration were plotted against the density ratio as 
shown in Figures 6-12 through 6-15. These velocities and accelerations 
were determined at a distance of 45 cm from the starting poing. The 
curves show that: 


1 . 


The lighter the sphere, the faster it will roll for a given slope. 
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VOLUME 
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SPHERE DENSITY, y ,®7 C e 

Fig. 6-10. Volume change as a function of sphere 
density, y , and soil density, y . 
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Fig. 6-12. Density ratio vs velocity and acceleration for spheres rolling on 
loose air-dry Yuma sand (cars 1 and 2). 
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Fig. 6-13. Density ratio /s velocity and acceleration for spheres rolling 
on dense air-cry Yuma sand (car 3). 




Fig. 6-14. Density ratio vs velocity and acceleration for spheres rolling 
on dense moist luma sand (cars 4 and 8). 




Fig. 6-15. Density ratio vs velocity and acceleration for spheres rolling on 
loose moist Yuma sa d (cars 5 and 7). 



2. The steeper the slope angle, the faster the sphere will roll 
for a given density ratio. 

3. The lighter the sphere and the steeper the slope angle, the 
greater will be the acceleration of the sphere. 

4. The acceleration of the sphere will be a relatively small 
fraction of the acceleration of gravity, and always less 
than the acceleration of the same sphere on a frictionless 
slope. 

From study of the curves of distance, velocity, and acceleration 
vs. time, an example of which was presented in Figure 6-11, it also 
appears that spheres tend to roll at constant acceleration down a slope 
provided soil and slope conditions are uniform. Of the some 200 spheres 
rolled during the experimentation reported herein, about 76% could be 
described as rolling at constent acceleration at 45 cm. after start. 

Velocity and Acceleration vs Sphere Diameter 

Velocity and acceleration as a function of sphere diameter are plotted 
in Figure 6-16 through 6-18. It appears that for the range of sphere sizes 
tested, velocity and acceleration are 

1. little effected by sphere diameter in dry sand, 

2. more effected by sphere diameter in dense moist sand, and 

3. distinctly effected by sphere diameter in loose moist sand. 

This implies that the dependence on sphere diameter may be related to 
cohesion which is nearly zero for dry sand, but has a small value for moist 
sand. 

However, it is important to note that the range of sphere sizes tested 
was very limited. Since lunar boulders may be about 10 meters in diameter, 
it may not be possible to extrapolate from the curves or Figures 6-16 
through 6-18 to such large sizes. If the curves of acceleration are 
extended to a 10 meter diameter sphere, the acceleration would be very 
small or zero. One may conclude, perhaps, that such large spheres would 
tend to roll at nearly constant velocity, at least much more so than small 
spheres. This implication could be significant since at constant velocity 



Fig. 6-16. Sphere diameter v > velocity and acceleration for spheres rolling 
on air-dry Yuma snd (cars 1, 2, and 3). 
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Fig. 6-17. Sphere diamete ■ vs velocity and acceleration for spheres rolling 
on dense moist Yuma sand (cars 4, 8, and 9). 





Fig. 6-18. Sphere diameter s velocity and acceleration for spheres rolling 
on loose moist Y ima sand (cars 5,6, and 7). 
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we have a unique relationship between slope angle, a, and the w/D 
ratio, as expressed by Equation (5-76). 


Slip vs Density Ratio, Sphere Diameter, and Distance Traveled 

Information on slip, defined as the distance revolved minus the 
distance traveled divided by the distance traveled, is presented in 
Figures 6-19 through 6-21. Figure 6-19 presents slip as a function of 
the density ratio. Figure 6-20 presents slip as a function of sphere 
diameter, and Figure 6-21 presents slip as a function of the distance 
traveled to the interval where the slip was determined. The following 
behavior is apparent: 

1. For loose moist sand, slip is primarily a function of the 
density ratio and little affected by the slope angle. 

2. For dense moist sand, slip is primarily a function of the slope 
angle and less affected by the density ratio. 

3. The data for slip in dry sand showed more scatter. However, 
the trend of the slope of the lines is similar to those for 
loose moist soil. 

4 J. xp dpptidl £> LU DC Cl tailiy UUUij-'J.CA j.uuwwxwa vj. W . 

For dense dry sand, slip decreased with increase in diameter. 
For loose moist sand (the other extreme tested) , slip increased 
with increase in diameter. 

5. Figure 6-21 shows that at least for travel distances greater 
than about 100 cm, slip was essentially independent of the 
distance traveled. This means that slip is essentially inde- 
pendent of sphere velocity and acceleration. This statement 
may also be true for a bouncing sphere? sphere E, car 8, 

a = 20° , bounced continuously beyond a travel distance of 
approximately 250 cm and yet the slip remained constant, as 
shown in Figure 6-21. 

6. All slip is negative for a freely rolling sphere? that is, the 
sphere travels further than its revolutions for rolling on a 
hard surface would indicate. 
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Fig. 6-19. Density ratio vs slip for spheres rolling on Yuma sand. 









DISTANCE REVOLVED - DISTANCE TRAVELED 
SL P * DISTANCE TRAVELED 


Fig. 6-21. Distance traveled vs slip for spheres rolling on Yuma sand. 
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7. Slip increases with the density of the sphere and the slope 
angle, as shown in Figures 6-19 and 6-20. 


Measurements of Shearing Soil Wedges 

Theoretical studies of the inertia forces associated with the 
shearing soil wedge showed that the soil acceleration is a function of 
the exponent, n. While the experimental studies were not aimed directly 
at the specific measurement of n, analysis of certain sequences of frames 
in the films permitted determination of the movement of a shearing soil 
wedge. The results of these measurements are presented in Table 6-4. 

Based on the data in Table 6-4 distance, velocity, and acceleration 
are plotted as functions of time in Figures 6-22 and 6-23. Velocities 
were determined from the slope of the distance curves, and accelerations 
were determined from the slope of the velocity curves. The experimental 
values of maximum velocity and acceleration as reported in Table 6-4 were 
determined from Figures 6-22. and 6-23. The experimental magnitude of 
sphere displacement (the distance through which the sphere moved in the 
same time) was determined from distance vs. time curves for the particular 
spheres . 

Theoretical soil wedge displacements were computed using the procedure 
suggested in Chapter 5 (x^^ 1 = z - z^) . The theoretical and experimental 
soil displacements reported in Table 6-4 are for comparable sphere dis- 
placements. Theoretical values for soil acceleration were computed from 

equation (5-60) using experimental values for v and t . The exponent. 

max so 

n, was assumed equal to one. 

If the data presented in Table 6-4 and Figures 6-22 and 6-23 can be 
considered representative of the movement of shearing soil wedges, certain 
inferences can be made: 

1. The exponent, n, appears to be close to one; i.e. the slopes of 
the increasing and decreasing sides of the velocity curve have 
nearly equal absolute magnitudes. 

2. The acceleration function assumed in Chapter 5 is undoubtedly 
an oversimplification. It appears that the velocity does not 
increase parabolically , as predicted by theory for n = 1, but 


Table 6-4. Measurement of Shearing Soil Wedges 


i H 

rH 

<d 

u 

II 

•H 

■P 

c 

0) 


H 

gj 

0 

6 

4) 

p 


to 

E-* 

to 

(d 


rH 


•H 

<d 

O 

a 

cn 

id 



II *H I 4H 
all rd 
C 

I 


•rl (d \ 

O G B 

cn id o 


J> 

in > o 


120.0 


576.0 





O 


CO 




00 


01 




• 


• 




rH 


O 



CN 





00 

00 





01 

• 





• 

o 





1— 1 

o 


o 


o 


• 


• 


• 


r- 


o 


o 


ID 


CN 


vD 




CO 


O 






rH 




rl ft *H 6 
-HO) O O 
O *H to 
m a X 


oi oo o o o o 
in O CO CO VD H 
O O CN O H 

• ••••• 
O O O O rH rH 


m^ioooo^rco 

oom^iocooHCM 

O O O O O rl rl rl 


o in o o 

00 rl 01 ID 

o cm oi r* co 
• • • • 
O O H M 



i 

* i 

a 

u 1 

•H 

<D 

Eh 

to 




>4* 

u J 

4 N 

i ' 

i — i 

V_J 

4. N 

>4* 


4JU 


IN 

>4' 

Vi-r 

IN 

^4’ 

vu 

U J 


cn 

o 

O 

O 

rH 

rH 

CN 

rH 

rH 

rH 

rH 

rH 

CN 

CN 

CN 

CN 

CO 

CO 

CO 

CO 

^i 

o 

rH 

rH 

t— l 

rH 

rH 

rH 

O 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

El 


co in id r-* co oi miDr-oooiOiHCNco 

CN CN CM CM CN CM CN rl rl H r| 


id r- oo oi o 

i— I rH rH i— I CN 



*With respect to mark on sand surfac 





















153 


I 

f 

K 


i 


s 

u 


lAj 

O 


tn 



l 


l 


o 

s* 

E 

« 


> 



UJ 

> 


o 



o 


o 

o 




Fig. 6-22. Soil movement as a function of time for 
sphere B, car 1 , a = 20°. 
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3. 

SUMMARY 
Testing 

The testing reported on in this chapter led to a consideration of 
certain specific relationships. Some of these relationships were investi- 
gated in order to increase information on rolling sphere — soil slope 
interaction; whereas , others were investigated in order to make a comparison 
with theoretical predictions. The relationships considered were: 

a. For general information 

1. Soil volume change associated with forming of th<» 
track, as a function of sphere and soil densities. 

2. Distance, velocity, and acceleration of the rolling 
spheres as a function of time. 

3. The relationship between velocity and the density ratio, 
and between acceleration and the density ratio for the 
rolling spheres. 

4. The relationship between velocity and sphere diameter, 
and between acceleration and sphere diameter. 

5. The relationship between slip and density ratio, sphere 
diameter, and distance traveled for the rolling spheres. 

b. For both comparison with theory and general information 

1. Slope angle for constant velocity rolling at minimum 
slope angle required for rolling. 


somewhat more slowly. As shown in Table 6-4, theoretical 

values of a were about twice the experimental values, 
max 

This can be explained from the difference in the theoretical 
and experimental velocity curves shown in Figure 6-22. The 
initial slope of the theoretical curves, is about twice 

the initial slope of the experimental curve, a'^ . 

However, the comparison does suggest that the theory gives 
results that are within a factor of two of the correct value. 


2 . 


The relationship between track width and density ratio 
for rolling spheres of constant size. 

3. The relationship between track width and sphere diameter 
for rolling spheres of approximately constant density. 

4. Distance, velocity, and acceleration of soil wedges 
moving during general shear. 


Comparisons with Theory 

For the items listed under b. above, the following conclusions can 
be made: 

1. The experimental data plotted in Figure 6-6 confirm the predicted 
relationship between track width and slope angle. 


= sin 2a, (5-76) 

for constant velocity rolling or for minimum slope angle required 
for rolling. This also confirms that the sphere — soil contact 
area must be approximately semicircular in a view normal to the 

<= 1 r\r>£5. 

jc — - 

2. The comparisons presented in Figure 6-7 confirm the validity 
of the proposed Equations (5-72) and (5-74) , 



0.188N + 

ys 



(5-72) 



for the range of sphere and soil conditions investigated. It 
is believed that the proposed equations are applicable to many 
other sphere and soil conditions as well. 

3. The comparison in Figure 6-8 is further evidence for the general 
validity of the proposed equations, particularly Equation (5-74) 


4. The data on soil wedges moving during general shear are 

inadequate for definite conclusions. The following inferences 
can, however, be made: 

1. The velocity does not increase parabolically for n = 1, 

but somewhat more slowly. It appears that for this 

reason theoretically predicted accelerations, a , may 

max 2 

be too high, 

2. The proposed theory will give answers for which are 
within a factor of two of the correct value. 

Overall Value of Testing and Theory 

In addition to the comparisons between theoretical predictions and 
experimental data, a considerable amount of experimental data on the 
dynamic aspects of rolling sphere — soil slope interaction were presented. 
This information is believed to be valuable for a better understanding 
of the problem. 

It is believed that the proposed theory is adequate for use as a 
remote reconnassance method for estimating soil conditions. For controlled 
testing and where q^ can be assumed to be negligable, the proposed theory 
and the rolling sphere relationship may give values for <p which are close 
to those determined from the results of triaxial tests. 
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Chapter 7. UTILIZATION OF DEVELOPED THEORY 

The implications of the theory, particularily the added inertia 
term in the bearing capacity equation, can now be considered. The 
relative contribution to total soil resistance of inertia, friction, 
cohesion, and surcharge are evaluated for selected situations in this 
chapter. 

A method of analysis of boulder tracks, based on the developed 
theory is outlined, and the lunar boulder tracks analyzed in Chapter 2 
sure reanalyzed by the proposed method. 

RELATIVE CONTRIBUTION OF INERTIA, FRICTION, COHESION, AND SURCHARGE TO 
TOTAL SOIL RESISTANCE 

The relative contribution of each of the terms in Equation (5-72) , 



0.188N 

ys 



+ 




(5-72) 


is evaluated below. 

The inertia term as expressed by Equation (5-70) is 


WY 

's 


b(z 


\ s 

s v; w 3 d' 


(5-70) 


where 


(1 + 2n) vi 

(sin-^y 9° 


B = 20.3 


(5-71) 
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The flow diagram below shows how the inertia term is computed: 


♦ 



2 ! 

wy 


Both z and can be determined from Figure 5-6. Knowing w/D, 

z would be determined by drawing a horizontal line from w v /D to z^/r as 

if these ratios were w/D and z/r. The value of z would be determined 

v 

U.. C d : 4-1. n At» /r\2 

voi. 4.11^ i. ' w n x k.n m m*wnu yuxuv. ^ ju u v / l/ • 

It was shown in Chapter 6 that the theoretical maximum acceleration 
may be about twice the experimental or actual value. Therefore, the 
values of B obtained from Figure 5-17 will be divided by 2. 

The other terms in Equation (5-72) are determined using Equations 

(5-20), (5-22), and (5-25), and Figures 2-4, 2-5, and 2-6. The reduced 

bearing capacity factors, N , N , and N , calculated from Equations 

ys c s qs 

(5-20) , (5-22) , and (5-25) would, in general require an increase in 

the friction angle of about 6 degrees from that required using , 

N , and N to get the same soil resistance, 
c q 

Since Equation (5-70) involves many parameters, all of which 

could vary independently in a real situation, any evaluation will be 

limited by the assumed conditions. It was considered desirable first 

to make an evaluation of the inertia term with variable sphere velocity, 

v ; and next to make an evaluation of the inertia term with variable 
s 

sphere diameter, D. 
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Consider a spherical rigid wheel 60 cm in diameter being pulled 
on a horizontal surface of sand having the following properties: 

2 2 

apparent cohesion = c = 1 kN/m ** 10.2 gm/cm 
friction angle = $ - 37.5° 

density = Y g = 1.6 gm/cc 

(Such values were reported by Costes and Mitchell, 1970, and Apollo 11 
Soil Mechanics Team, Costes et al., 1970, from their examination of 
Apollo 11 samples.) If the wheel load is varied leading to varying 
penetrations or w/D ratios and the velocity is varied, the magnitude 
of the inertia term will vary as shown in Figure 7-1. 

\ 

Figure 7-la shows the percent contribution of the inertia term, 

q /wy , to total soil resistance as a function of velocity. For these 
x s 

lunar soil conditions inertia contribution begins to be important at 
a velocity of about 5 miles per hour. For soil conditions used in 
testing the instrumented spherical wheel on Yuma sand at a velocity of 
about 1.3 mi/hr or 60 cm/sec, the resistance due to soil inertia was 
determined to be 3.4% — relatively insignificant. 

Figure 7-lb shows the relative contribution to total soil 
resistance from each of the terms in Equation (5-72) at a selected 
sobers velocity of 16 km/hr. For other wolnri>iPB. t-.he magnitude of 
the density, cohesion, and surcharge terms would remain unchanged but 
the inertia term would vary as shown in Figure 7-la. Consequently, a 
change in the percentage contribution would result. The striking 
feature of Figure 7-lb is the contribution due to cohesion. For small 
loaded areas, w/D, the cohesion term provides most of the resistance 
in spite of the fact that the assumed value of unit cohesion was very 
small. The reason for this is that a large value of 4> gives a high N^. 

Next we are to consider an evaluation of the inertia term with the 
size of the sphere as the principal variable. This evaluation is much 
more complex because as the size of the sphere increases, penetration 
and w/D ratio increase as shown in Figure 6—8 . As the w/D ratio 
increases the slope angle required for rolling increases. Because of 
these dependencies, the velocity of the sphere is, in fact, a function 
of the size of the sphere. Some experimental data are available for 
small spheres (7 to 25 cm in diameter) . However, no data are available 


TOTAL 


a) Resistance due to soil inertia as a function of sphere velocity 



b) Contribution to total soil resistance as a function of w/D 

Fig. 7-1. Relative contribution to total soil resistance from 
friction and density, cohesion, surcharge, and soil 
inertia for a spherical wheel ( D = 60 cm). 
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on the velocity of large (10 meter) rolling lunar boulders. It is also 
doubtful that direct experimental verification will ever become 
available for large boulders since it is impractical to conduct tests 
on 10-meter spheres. Therefore, in determining B, it is at the present 
time difficult to assess values for the velocity term, v 2 /gD, for large 
boulders . 

We begin by plotting data for smaller spheres and extrapolate to 
larger spheres. The velocity term was determined for a 20° slope from 
Figures 6-16 and 6-17 for dense dry Yuma sand (car 3) and dense moist 
Yuma sand (car 4) . These Figures indicate that larger spheres roll 
more slowly them smaller spheres, and that a 10-meter sphere rolling 
on the same soil would roll very slowly indeed. The values of v 2 /gD 
obtained are plotted as the experimental points in Figure 7-2. 

The track width for many of the lunar boulders reported in 
Chapter 2 was determined at a location a short distance before the point 
where the boulder came to rest or in a region where it appeared to have 
been rolling uniformly. A velocity of about 2.5 meters per second may 
be representative for these boulders. This value will be assumed . 

This would give v 2 /gD = 0.08. Using this value as one of the points 
in Figure 7-2, a curve was fit to the data and is shown in Fiqure 7-2. 

For selected values of D, the velocity term, v 2 /gD, was obtained 

from Figure 7-2; and the contribution of each of the terms in 

Equation (5-72) was determined. The same soil conditions (c = 1 kN/m 2 , 

4> = 37.5°, and y =1.6 gm/cc) were again used. The results are shown 
s 

in Figure 7-3. Figure 7-3a shows the relative contribution of each 
term in Equation (5-72) in earth gravity. In lunar gravity, the 
inertia and cohesion terms remain unchanged while the density and 
surcharge terms are divided by six. Figure 7-3b shows the relative 
contribution of each term in lunar gravity. 

As shown in Figure 7-3, the contribution to total soil resistance 
from cohesion is much influenced by the size of the sphere. For the 
small spheres (7-25 cm in diameter) tested at the WES, the moist Yuma 
sand behaved much as a cohesive soil, while the air dry Yuma sand 
behaved much as a cohesionless soil. That the theory predicted 
correctly the track widths for both conditions, as shown in Figure 6-7 
supports its general validity. 
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a) Earth gravity. 
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b) Lunar gravity. 


Fig. 7-3. Relative contribution to total soil resistance from 
cohesion, surcharge, friction and deputy, and soil 
inertia as a function of sphere diarrjter. 
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Figure 7-3b shows that for the conditions assumed, soil inertia may 
contribute about 20% of the total soil resistance during lunar boulder 
track formation. It will be recalled that the 20% was calculated using 
an assumed velocity of about 2.5 meters per second. Some lunar boulders 
probably rolled faster and some rolled more slowly. The important 
conclusion is that where as on earth the inertia contribution would have 
been relatively insignificant, on the moon the inertia contribution is 
significant. 

METHOD FOR SOLVING FOR $ FROM LUNAR BOULDER TRACKS 

The flow diagram on the following page shows how <p can be 
determined. 

For selected values of w, w/D, and a, and lunar soil conditions 
(c = 1 kN/m 2 , Y = 1.6 gm/cc, and y /y = 2) , the friction angle was 
determined as proposed. The resulting curves are presented in 
Figures 7-4 and 7-5. These figures can be used to solve for the 
friction angle of lunar soil and to study variability of lunar surface 
soils. 

In evaluating lunar soils by the proposed method, certain 
behavior accompanying compressibility should be recognized. From 
studies with simulated lunar soil, Houston, Namiq, and Mitchell (1970) 
showed that the density , cohesion and friction angle increase with 
depth. When a friction angle is determined from the boulder-track 
relationship, it is an average value for the particular failure surface. 
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REANALYSIS OF LUNAR BOULDER TRACKS BY PROPOSED METHOD 
Analysis and Results 

For the boulder tracks analyzed in Chapter 2, the assumed soil 
parameters were: 

c = 0.1 kN/m 2 = 1.02 gm/cm 2 

y = 1.6 gm/cc 
s 



The approximate average of the measurements of the w/D ratio from 
the 69 boulder tracks was about 0.69. The average slope angle was 
about 15 degrees, and the average track width was about 6.3 meters. 

On the basis of all the boulders analyzed , the average friction angle, 
4>, worked out to be 34 degrees. 

Tests performed on the lunar samples brought back to earth with 
Apollo 11 permit a refinement of the soil parameter values. The 
following values are more reasonable: 

c = 1.0 kN/m 2 = 10.2 gm/cm 2 

Y s = 1.6 gm/cc 



With these soil and rock parameters, the 69 lunar boulder tracks were re- 
analyzed using the proposed theory (Figures 7-4 and 7-5) . The results 
are presented in Table 7-1 and Figure 7-6. As shown in Figure (7-6) , 
most of the values of <}) were between 24 and 47 degrees. The values 
ranged from 19 to 53 degrees, and the arithmetic average of all values 
of <J) worked out to be 37 degrees. 

The values of <}> were also plotted against the boulder diameter 
as shown in Figure 7-7. This figure includes all results for all 
w/D ratios and slope angles a. It can be concluded that the higher 
values of <}> were associated with the larger boulders. This seems to 
imply that 

1. The friction angle for the lunar regolith increases with 


Table 7-1A. Results of boulder track analysis, Orbiter II photographs. 
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Table 7-1B. Results of boulder track analysis, Orbiter III photographs. 
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Table 7-lc. Results of boulder track analysis, Orbiter v photographs. 
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Fig. 7 - 7 . Results of <j> plotted against boulder diameter. 
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depth. This is probably a consequence of an increase in 
soil density with depth. These inferences are based on the 
fact that the larger boulders would generate a deeper failure 
surface. Analogous behavior has been predicted from studies 
on simulated lunar soil (Houston, Namiq, and Mitchell, 1970). 

2. Alternatively, it is possible that the larger boulders were 
influenced by a layer of harder material below the lunar , 

regolith. This could have resulted in a lower value of w/D 
and an apparently higher «J>. In this case the indication is 
that the thickness of the lunar regolith in the vicinity of 
the boulder tracks was close to the radius of the larger 
boulders (10 to 15 meters) . 


Discussion 

As shown in Table 7—1 and on Figure 7-6, <p varied considerably - 
from 19 to 53 degrees. In some locations, the measurements indicated 
a significant variability within the same slope. That is, the sinkage, 
z/w, of one boulder was significantly greater than that of an adjacent 
boulder also analyzed (Frames III-125H, V-95H, and V-168H) . Possible 
interpretations for such variations would be that : , l) the soil properties 

. . ... _ V 

airrffrfln gx T!!*? TWO JLOC3.Tll.OnS , tltlO JDOUJLCL92T C3»l2Singf tilS w2T3.Ck SOITiChCV 

did not conform to the assumptions of the theory used; perhaps it was 
not spherical or it had inertial energy beyond that assumed, and 3 ) 
other as yet unknown causes. For most of the locations, on the other 
hand, the results were more constant. 


For many tracks, the measured track width, w, was equal to the 
boulder diameter, D. Hence, we can conclude that in such cases the 
thickness of unconsolidated soil on the slope must have been at least 
half the boulder diameter. 

Since the Surveyor results (Surveyor VII Mission Report, 1968) 
and Apollo 11 results (Costes and Mitchell, 1970) quite consistently 
indicate that the friction angle for lunar soil may be in the range 
of 37 to 39 degrees, it is of interest to consider possible causes 
of the considerably larger range found here. As stated before, the 
results of this analysis indicate an average 4> of 37 degrees with a 


range from 19 to 53 degrees. It is appropriate, however, to assess 
the degree of confidence that may be placed in these results. 


Sensitivity of results to variations in soil density: 

The influence of variations in densities of rock and soil can 
be appreciated by considering Equation (5-75) , which gives the left 
side term of Equation (5-73) . The q m /wy s term is directly proportional 
to the density ratio, y /y , and inversely proportional to the w/D ratio 

IT S 

cubed. As shown in Figure 2-7, for any one w/D ratio, the term 

could almost triple within the range of possible values of y /y between 

XT S 

1.7 and 3. This range of values for the q^/w Y g term would lead to 
a ± 5° change in the apparent friction angle, <j>, determined from the 
boulder track relationship. 


Sensitivity of results to variations in cohesion: 


The significance of cohesion depends on the relative magnitude of 

the cohesion term in Equation (5-73) . The sensitivity of the results 

to variations in cohesion was investigated using average values (<J> = 37° , 

c = 10.2 gm/cm 2 , z/w = 0.234, w = 7m, y = 1.6 gm/cm 2 , and a = 15°) 

s 


for the boulder tracks analysed. For a 10 and 100 told increase in 
cohesion, the required friction angle for the same soil resistance is: 


Cohesion 
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c 


<J> 


H* 

O 

u> 

dynes/cm 2 

39 

10 4 

dynes/cm 2 

37 

10 5 

dynes/cm 2 

27 


Sensitivity of results to measurements: 

To see how serious errors or inconsistencies in measuring the 
boulder and track dimensions might be, the measurements were made on 
two different scale photographs as previously described. The difference 
in the two measurements were computed as follows: 
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For the diameter 
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For the track width 
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D larger 
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smaller 

w 

larger 
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X 100 


x 100 


Similarily, the difference in the w/D ratio was computed by 


A 



(w/D) 

'(w/D) 


smaller 

larger 


x 100 


The results of the comparison and the resulting change in the 

friction angle, A<J>, are shown in Table 7-2. The basis for this 

comparison is (w/D) = 0.69 and d> = 37°. 

ave ave 


Table 7-2. Sensitivity of results to measurements 


Comparison 

a*? 
< — 

Aw 

{%) 

A w/D 
(%> 

A<)> 

(°) 

Average 

9 

12 

9.6 

2 

Maximum 

32 

28 

22 

5 


As shown in Table 7-2 above, the average difference in <f> for the two 
measurements would be about 2 degrees. If it can be assumed that the 
friction angle resulting from averaging the measurements, as was done 
for the results presented, is closer to the correct value than the 
result of either of the measurements independently, the friction angles 



presented in Table 7-2 may be regarded as including a one degree 
average error and about 3 degree maximum error due to inconsistencies 
in measurements. This is quite good considering the scale of the 
photographs; as better photographs become available, it should be 
possible to considerably reduce measurement inaccuracies. 


Sensitivity of results to the slope angle: 

The influence of changes in the slope angle, a, on the resulting 
friction angle was determined from Figures 7-4 and 7-5. Again the 
basis for the evaluation was the average friction angle of 37 degrees, 
and the average slope angle of 15 degrees. The results are presented 
below: 


Slope angle, a 

0 ° 

5° 

10 ° 

15° 

20 ° 

25° 

*5hO 


Friction angle, <J> 

33° 

34° 

35.5° 

37° 

38.5° 

40° 


Since slope angles cannot be very accurately determined at this time, 
the results of <p presented in Table 7-1 should be considered to have 
a probable error of ±1.5° due to incorrect estimation of the slope 
angle. This would correspond approximately to a ±5° difference in 
the slope angle as shown above. 


Summary 

From a consideration of the possible effects of incorrect 
assumptions as to soil parameters, incorrect measurements, incorrect 
estimates of the slop* angle, and inadequacies of the theory, probable 
uncertainties in the results are as follows: 


Cause : 


Effect on (j) (°) 


Variation in density ratio from 
Y r /Y g = 1 to 3 

Variation in c (for every 

1 x 10 4 dynes/cm 2 , 10 3 < c < 10 5 dynes/cm 2 ) 

Incorrect value of w/D ratio 
(measurements ) 

Incorrect slope angle (±5°) 


±5 

±1 to 2 
±1 to 2 
±1 to 2 


Of the causes listed above, incorrect estimates of the density 
ratio have the largest effect on the results. Since the values of 
<P varied from 19 to 53 degrees, it appears that the density ratios 
may have varied considerably for the boulder tracks analyzed. Possibly 
the state of soil compaction or the specific gravities varied. Incorrect 
assessment of cohesion and determination of the w/D ratio and the slope 
angle probably also contributed to the spread in the values of <j>. 

CONCLUSIONS 


Certain implications and applications of the developed theorv were* 
considered in this chapter. The following conclusions can be made: 

1. The cohesion term provides most of the resistance for small 

(5 x 5 cm to 10 x io cm) loaded areas in moist sand; because 

although c may be as low as 10 gm/cm 2 , N is large for a 

c 

high friction angle. 

2. The resistance due to inertia of the moving soil is more 
significant on the moon than on earth. 

3. A method for solving for the friction angle for lunar soils 
from the boulder track relationship was outlined, it is 
believed that Figures 7-4 and 7-5 can be used to determine 
a friction angle and to study variability of lunar surface 
soils. Provided assumptions regarding other soil and rock 
parameters are realistic and measurements of boulder and 
track dimensions are adequate, the proposed method should 
give friction angles comparable to (f> determined from 
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conventional triaxial tests. 

4. The 69 lunar boulder tracks reviewed in Chapter 2 were 
reevaluated using the proposed method. The values of <p 
ranged from 19 to 53 degrees, and the arithmetic average for 
the friction angle of lunar soil was found to be 37 degrees. 

5. Lunar soil and rock properties (cohesion, density, and 
friction angle) are possibly more variable than anticipated 
from previous investigations. 

6. The results also indicate that the soil parameters (density, 
cohesion, and friction angle) increase with depth below the 
lunar surface (Figure 7-7) . 
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Chapter 8. CONCLUSIONS AND RECOMMENDATIONS 


GENERAL 

The studies reported herein were motivated by observations of 
boulder tracks both on the moon and on earth. From such boulder-track 
records usually only the track width, the boulder diameter, and the 
slope angle can be measured. With this in mind, both experimental and 
theoretical studies were aimed at relating soil and rock parameters 
to the measurable w/D ratio. 

Provided other soil and rock parameters are known from independent 
analysis, it is possible to determine one unknown from the boulder-track 
relationship. The experimental data and comparisons presented in 
Chapter 6 show close agreement with theoretical predictions. The 
friction angle of lunar soil, for example, was found to be 37 degrees 
using the proposed theory. Independent studies of lunar soil from 
Surveyor tests and Apollo 11 samples have resulted in nearly the same 
value. It is believed that soil conditions can be evaluated from the 
boulder-track relationship. 

Xn 2T3. , tiis s tudis s provided v^lusbls i-nt 01 rm 3 .tj.cn on ttc 
failure mechanism associated with rolling sphere-soil slope 
interaction. The failure mechanism for wheel-soil interaction is 
believed to be similar at least for nearly rigid wheels. The 
theoretical developments are a first attempt to incorporate the 
influence of soil inertia as a part of soil resistance. 

THE FAILURE MECHANISM 

As a sphere rolls down a sand slope, the most important 
characteristics of behavior are: 

1. Track formation is accompanied by soil volume changes and 
general shear. 

2. For loose sands and at the beginning of rolling in general, 
soil deformations consist of volume change. 

When volume change (compression) can no longer account for 


3. 
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4. 


5. 


6 . 


7 


8 . 


9. 


all the soil that must be displaced, shear planes develop. 
Within the shear zone in front of the sphere there may be 
about two shear planes per cm of depth of the shear zone. 

The shearing is a discontinuous process although rolling 
appears to be continuous to the naked eye. Separately 
identifiable shearing soil wedges are pushed out in front of 
the sphere. These failure wedges appear to form where the 
soil is weakest; sometimes directly in front of and sometimes 
slightly to the sides of the sphere. 

The maximum pressure beneath the sphere is about twice the 
average pressure, and the maximum pressure appears to control 
the development of shear failure. 

The contact pressure distribution is bell-shaped and can be 
closely approximated by a parabolic function. 

Vertical sections within the soil are bent forward. This 
forward movement is pronounced; very little lateral 
movement takes place. 

The magnitude of the resultant of soil reaction can usually 
be considered equal to the weight of the sphere. The 
direction of the resultant can be assumed vertical for 
prdCuiCal purposes • 


THEORY 


Equations (5-70) through (5-76) can be used to solve boulder-track 
problems. These equations are quite complete in terms of consideration 
of the many variables involved in the rolling sphere (boulder) -soil 
slope interaction problem. Because of the number of variables, 
however, many assumptions were made. Further, due to the assumptions 
and the nature of the theoretical development, most of the equations 
are semiempirical . Independent checks of these assumptions by 
experimentation are recommended. However, it is believed that the 
proposed equations can provide a valuable basis for a remote 
reconnaissance technique for estimating approximate values of soil 
parameters . Provided assumptions regarding other soil and rock 
parameters are realistic and measurements of boulder and track 
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dimensions are adequate, the proposed method should give friction angles 
comparable to <}> determined from conventional triaxial tests. 


SOIL PROPERTY EVALUATION 


From Rolling Spheres 

Based on the investigations reported herin, it appears that it 
should be possible to study surface soil conditions on the basis of 
sphere or spherical wheel-soil interaction. It has been noted that 
the resulting track width is very sensitive to changes in soil 
conditions. Soft pockets are easily detectable from a sudden increase 
in track width. Data obtained by pulling a spherical wheel on the 
soil surface would be required. By monitoring the velocity and 
measuring or continuously photographing the resulting track, all the 
parameters in Equation (5-72) could be systematically studied. 


From Rolling Boulders 

Velocity information is not directly available from boulder-track 
records. Further, assumptions must be made as to other soil and rock 
properties such as density and cohesion. Therefore, the answers 
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assumptions. Nonetheless, it is believed that the study of boulder- 
track records can provide valuable information on the variability of 
soil conditions and estimates of soil parameters for unknown areas. 


Properties of the Lunar Regolith 

The lunar boulder tracks were reevaluated using the method 
proposed in Chapter 7: this method is based on the improved theory 
developed in Chapter 5. The following conclusions seem to be appropriate 
from the results of the 69 lunar boulder tracks investigated: 

1. The results of <j> ranged from 19 to 53 degrees with the 
majority being between 24 and 47 degrees. 

2. The average friction angle, <f>, was found to be 37 degrees. 

3. Based on the range of 4>, lunar soil conditions appear to be 
quite variable. 

, Soil density, cohesion, and friction angle probably increase 


4 


with depth below the lunar surface. 


Possible implications of the range of <f> reported above for the 
boulder tracks from different areas of the moon are: 

1. The state of soil compaction of lunar soil may vary. 

2. The cohesion may vary. 

3. The friction angle may vary. 

4. The measured track width, whiab is the basis for the range 
of <j>, was not representative of the crest-to-crest track 
width immediately after the track was formed. 

Of the above possibilities, undoubtedly some of the track widths 
were not representative of the crest-to-crest track width immediately 
after the track was formed. On the other hand, some variability of 
soil conditions should be expected. Particularly the state of 
compaction of the lunar regolith may vary. 

In applying the proposed theory to the solution of lunar boulder- 
track problems, it should be remembered that inertia forces are 
relatively more important on the moon than on earth. For the 69 lunar 
boulder tracks described in this report, soil inertia was considered 
provide 20% of the total soil resistance. 

RECOMMENDATIONS 

It is recommended that the boulder track method be considered as 
a remote reconnaissance technique for preliminary study of soil 
conditions at different points on the moon and on other planets 
where bouldee tracks are observed. 

It is recommended that efforts in the near future be directed 
also to a study of wheel-soil interaction. It is believed that the 
failure mechanism described is applicable to many wheel-soil 
interaction problems. With some modifications the proposed theory may 
be suitable for the study of wheels. Since much data are already 
available in the literature on wheel-soil interactions, it might 
be possible to check the proposed theory without extensive testing. 
Such studies are recommended. 
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The failure mechanism proposed should be considered in relation 
to existing trafficability theory. A spherical wheel could be designed 
and tested to study the inertia effects. The testing of such a wheel 
could provide information on the exponent, n f and on the acceleration 
of the shearing soil. This could be of considerable value in the 
design and performance prediction of terrain vehicles and lunar rovers. 
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semicircular contact area based on crest-to-crest track 
width , w 

longitudinal section area of failure wedge for a footing 

longitudinal section area of failure wedge for a sphere 

total cross-sectional area of moving soil including 

cross-sectional area of track below crest of track 

linear acceleration of sphere or boulder 

final acceleration (deceleration) of shearing soil 

initial acceleration of shearing soil 

maximum acceleration of shearing soil 

acceleration of shearing soil 

abreviation used in solving Bquation (5-61) 

integration constants 

width of equivalent rectangular area 
approximate width of shearing soil wedge 
apparent cohesion 

apparent cohesion determined from in-situ tests 

diameter of sphere or boulder 

relative density of sand in percent 

drawbar puli (towing force) 

drawbar pull corrected for carriage inertia 

dimension defining triangle base proportional to d Q 

dimension of cone underneath sphere 

resultant of soil reaction 

soil inertia force 

towing force for sphere 

component of F parallel to slope 

component of F normal to slope 

cone index gradient (cone penetration resistance gradient) 
acceleration of gravity 
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moment of inertia of sphere = - 'mr 
kilo newtons per meter squared 
coefficient of subgrade reaction for sand 
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length of equivalent rectangular area 
mega newtons per meter cubed 
mass of sphere or boulder = W/g 
mass of shearing soil wedge 
newtons 

Meyerhof (1951) bearing capacity factors 
bearing capacity factors for rolling spheres 
exponent defining acceleration function (n = - a ^/ a f) 
pressure cell 

unit bearing capacity, bearing pressure, or average soil 
pressure 

unit bearing capacity in earth gravity 

unit soil resistance due to inertia of moving soil 

unit bearing capacity in lunar gravity 

unit soil resistance due to general shear 

unit soil resistance due to volume change 

surcharge pressure 

resultant of wheel load and towing force 
radius of sphere or boulder 

distance from center of sphere to line of action of P 
area of failure auricle for a footing par unit depth 
perpendicular to the paper 

area of failure surface for a rolling sphere per unit 
depth perpendicular to the paper 

shape factors (for shape of loaded area) for the bearing 

capacity equation 

time 

time of soil movement 
angular- acceleration of sphere 

volume of track rims above original soil surface, per 
unit depth 

volume of shearing soil wedge 
total volume ( | | + | v w |) 

volume of center of track below original soil surface , per 
unit depth 

velocity of sphere parallel to slope 


v 



w 

WL 


w 
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average soil velocity w.r.t. time 
maximum soil velocity 

velocity of soil moving along failure surface 
weight of sphere or boulder 
wheel load 

crest-to-crest track width 

crest-to-crest track width for maximum volume change 

before general shear (w^ < w) 

distance traveled by sphere (Chapter 6) 

coordinate directions 

soil displacement during general shear 

total displacement of a shearing soil wedge 

coordinate distance from location of maximum pressure 

perpendicular to direction of maximum pressure 

sinkage or track depth based on crest-to-crest track width 

track depth based on or track depth for maximum 

volume change before general shear 
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slope angle 

angle defininq equivalent free surface on Meyerhof's charts 
angle between vertical and direction of resultant of soil 
reaction for a sphere 

volume of compressed soil per unit depth perpendicular to 

track (AV = V - V ) 
w r 

maximum deflection of diaphragm 

maximum strain in diaphragm 

* dry density of soil in earth gravity 

soil density in lunar gravity (y = y /G) 

m s 

rock or sphere density in earth gravity 
soil density in earth gravity 

angle at sphere center defined by bisecting the track 
width, w 

angle at sphere center defined by bisecting w v 

angle of intersection of slope and upper end of passive 

Rankine zone 

angle of intersection of slope and failure surface 
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♦t 
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II 

III 


pressure at the sphere-soil contact parallel to the 
resultant, R 

maximum contact pressure at y = 0 
normal stress 
major principal stress 
minor principal stress 

angle of shearing resistance or internal friction 
angle of internal friction based on in-situ tests 
secant friction angle 
tangent friction angle 

base angle (angle defining frictional resistance at a 
footing base-soil contact) 
shear stress 

moisture content (percent moisture per dry-weight of soil) 
Orbiter II 
Orbiter III 
Orbiter V 


V 
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Appendix A. PREPARATION OF SOIL CARS AT THE 
WATERWAYS EXPERIMENT STATION 

Preparation of a soil car began by raking the soil as shown in 
Figure A-l to break up the previous soil structure. 

Next the surface was leveled with the top of the car by screeding 
with a straight aluminum member. Then the soil was compacted to desired 
density using vibrators attached to a heavy piece of plywood, as shown 
in Figure A-2. The WES personnel have developed techniques through 
experience by which a desired density can be produced and reproduced. 

Desired densities in dry (air-dry) sand can be achieved quite easily. 

Moist sand requires greater precision and care. 

The soil strength or resistance to penetration is tested by the cone 
index test as shown in Figure A-3. A photoelectric eye located at the 
top of the cone indicates contact with the soil and starts the automatic 
recording process. A load cell measures soil resistance. The slope of 
the stress vs. depth plot is the cone index gradient G expressed in MN/m 3 
in this report. In air-dry sand, the density would be determined from 
the curves in Appendix B, entering the charts with G determined as 
explained above. 

For moist sand, the density must be determined by an in-place density 
test. An in-place density box designed by Dr. Hvorslev was used for this 
purpose as shown in Figure A-4. Entering the charts in Appendix B with G, 
the density, and moisture content (determined as explained above) other 
parameters, such as cohesion, can be obtained. 

In general, the raking and compacting process produced a relatively 
uniform soil to a depth of about 20 cm. The failure planes for all of the 
spheres are believed to have been much shallower, and therefore, the 
preparation of the soil cars was probably adequate. That is, the underlying 
denser sand probably had little influence, on the results. 
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Fig. A-l. Raking the soil. 



Fig. A-2. 


Vibrating the soil . 
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Appendix B. PROPERTIES OF YUMA SAND 

(Determined by the U. S. Army Engineer Waterways Experiment 
Station, Vicksburg, Mississippi) 

Soil property relationships for Yuma sand as determined by WES are 
presented in Figures B-l through B-10. 

As shown in Figure B-10 there is generally a seven to eight degree 
difference between the angles of internal friction determined from 
laboratory tests and those determined from plate in situ shear tests. 

This difference is probably due to the different failure mechanisms 
associated with the various tests (Freitag, Green, and Melzer, 1970) . 

The main component of the plate in situ shear test device was a 
hollow aluminum plate 60 cm long, 30 cm wide, and 2.5 cm thick. Sand 
was glued to the surface to be in contact with the soil to assure soil-to- 
soil shearing. The plate was loaded vertically and pulled horizontally 
at a constant gpeed of 1 mm/min until failure occurred. Further details 
of the testing procedures for the relationships presented in Figures B-l 
through B-10 are to be found in the report by Freitag, Green, and Melzer 
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Fig. B-6. Results of vacuum tri axial tests 



















Fig. B-9. Results of in situ shear tests on air-dry sand 






Appendix C. PRESSURE CELL DATA 


This appendix contains the pressure cell data obtained from testing 
the instrumented spherical wheel at WES. The design of the wheel and 
the gages, as well as the testing procedure, was described in Chapter 4. 

As the wheel rolled through the soil , pressures registered by the 
gages were automatically recorded on the oscillographs presented in 
Figures C-l through C-5. These oscillographs provided the basic data 
for the pressure contour diagrams presented in Figures 4-7, 4-9, 4-11, 
4-13, and 4-15. 

In addition, the signals from the gages were fed into a computer, 
which selected certain points and filtered out the 60-cycle disturbances 
(small wiggles) shown on the oscillographs. These data are presented in 
Tables C-l through C-5. 
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Fig. C-2. Oscillograph, wheel load = 1241 newtons. 
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Computer-Selected Filtered Data 
Table C-l 

Wheel Load = 635 newtons 


Approximate Averages for Controlled Parameters 
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Explanation: 


Location as identified on lower edge of oscillograph 
Db— roa * Drawbar pull corrected for carriage inertia 

N = newtons 2 2 

All pressures are in kilo newtons meter (kN/m ) . 
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Table C-5. (Cont'd) 
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Table D-l. Track data 
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